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Abstract 

Finite  element  analysis  has  been  used  as  a  design  tool  for  many 
years,  with  structural  reliability  being  ensured  through  use  of 
a  liberal  factor  of  safety.  Unfortunately,  the  safety  factor 
is  a  blanket  insurance  against  all  hazards,  and  a  designer  has 
no  way  to  optimize  a  structure  against  any  particular  hazard. 
This  is  particularly  troublesome  in  the  fields  of  aero/astro 
design,  where  every  bit  of  mass  must  serve  to  maximum  utility. 

The  method  of  Stochastic  Finite  Element  Analysis  allows  a 
designer  to  model  any  loading  or  hazard  condition  as  closely  to 
reality  as  desired  by  using  an  appropriate  probability  distri¬ 
bution  function.  Through  a  Monte  Carlo  simulation,  the  finite 
element  model  is  subjected  to  the  probability  functions.  The 
cumulative  output  is  analyzed  for  trends  in  failure  probability 
and  the  design  is  altered  to  enhance  its  reliability,  repeating 
the  process  until  the  desired  level  of  reliability  is  achieved. 
The  resulting  design  is  optimal  for  the  imposed  conditions,  and 
compared  to  a  structure  designed  with  a  traditional  factor  of 
safety  approach,  is  either  lighter  or  more  reliable.  This 
demonstration  revealed  that  for  similar  reliabilities,  a 
stochastically  designed  wing  was  20%  lighter  than  a  wing 
strengthened  by  the  factor  of  safety. 

The  major  drawback  in  applying  the  method  of  stochastic  finite 
element  analysis  is  that  very  large,  complex  models  can  require 
extraordinary  amounts  of  computer  resources.  '  ■  . 


A  DEMONSTRATION  OF  THE  METHOD 
STOCHASTIC  FINITE  ELEMENT  ANALYSIS 


Introduction 


Traditional  structural  design  has  evolved  constantly  over  the 
years.  For  example,  early  cathedrals  were  designed  much 
stronger  than  required  for  structural  soundness  because  their 
designers  didn't  understand  structural  mechanics  and  used 
extreme  overdesign  to  avoid  any  possibility  of  structural  col¬ 
lapse.  (5:259)  As  engineers  came  to  understand  the  mechanics 
of  structures,  they  were  able  to  reduce  the  size  of  supporting 
structure  without  reducing  its  ability  to  safely  carry  a  load. 
Nevertheless,  because  of  uncertainties  in  material  properties 
and  loading  conditions,  as  well  as  inexact  methods  of  modeling 
reality,  structures  were  always  designed  to  be  stronger  than 
needed.  This  extra  capacity,  or  Factor  of  Safety  (FS)  is 
routinely  applied  in  structural  design  as  insurance  against 
failure  due  to  both  predictable  and  unpredictable  loads  and 
structural  characteristics.  Nearly  every  civil  work  cons¬ 
tructed  in  the  United  States  is  built  according  to  a  legislated 
safety  factor,  as  dictated  by  such  codes  as  the  Uniform 
Building  Code  (7:455)  and  the  Manual  of  Steel  Construction 
(2:Sec  2,21). 

Use  of  the  safety  factor  has  worked  well  primarily  because  it 
is  simple  to  apply  and  it  works  well  in  providing  a  sound. 
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stable  structure.  It  provides  a  reasonable  means  to  account 
for  the  perhaps  dozens  of  variables  that  might  affect  the 
ability  of  any  given  structure  to  withstand  a  required  load. 
The  designer  has  only  to  reduce  the  maximum  allowable  stress 
of  the  structural  material  by  the  factor  of  safety  and  pro¬ 
ceed  with  the  design  and  analysis. 

Is  there  any  reason  to  abandon  this  simple,  reliable  means  of 
ensuring  satisfactory  performance?  If  weight  is  a  concern, 
then  the  answer  is  yes.  In  the  quest  for  higher  performance 
aircraft  and  greater  useful  orbital  payloads,  every  pound 
within  the  load  bearing  structure  should  ideally  be  stressed 
to  just  short  of  failure.  Any  lightly  stressed  structural 
element  contains  excess  structural  weight  and  does  not  repre¬ 
sent  optimal  design. 

Obviously,  through  careful  design  and  the  use  of  a  factor  of 
safety  of  one,  this  result  can  be  obtained.  A  problem 
occurs,  however,  when  there  ar-  undetected  structural  de¬ 
fects.  Perhaps,  the  material  is  just  a  little  weaker  than 
planned  for,  or  the  load  profile  might  be  slightly  more 
active  than  anticipated.  A  design  engineer  should  be  able  to 
compensate  for  these  potential  failures,  however  rare  they 
might  be,  without  resorting  to  the  broad  brush  approach  of 
the  factor  of  safety?  Additionally,  suppose  the  designer 
assumes  that  some  part  of  the  structure  contains  a  flaw  or 
that  a  localized  structural  failure  has  occured  and  wants  to 
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assess  its  affect  on  the  rest  of  the  structure.  Obviously, 
then,  an  engineer  must  account  for  these  hidden  weaknesses, 
design  uncertainties,  and  assumed  damages  with  an  optimal 
design  and  still  obtain  a  robust,  durable  structure. 

If  he  can  estimate  what  is  the  chance  of  the  occurrence  of 
any  given  event,  then  he  can  plan  his  design  accordingly.  If 
he  knows  that  only  one  weld  of  every  ten  thousand  is  defec¬ 
tive,  and  the  structure  has  only  ten  welds,  he  can  design 
appropriately.  Concurrently,  if  he  knows  that  the  occurence 
of  even  one  bad  weld  will  precipitate  complete  structural 
collapse  and  loss  of  life,  he  can  adjust  the  design  to  com¬ 
pensate  for  that  too.  It  is  this  situation  which  brings  us 
to  the  foundation  of  this  thesis. 

A  body  of  statistical  mathematics  along  with  standard  finite 
element  analysis  techniques  have  been  combined  into  a  method 
called  stochastic  finite  element  analysis.  The  method  is 
completely  computer-oriented  and  relies  on  the  computer's 
ability  to  do  lengthy  finite  element  calculations  quickly  and 
repetitively.  The  recent  advent  of  cheap,  high-speed 
computer  time  makes  this  technique  feasible  because  of  the 
method's  potential  for  using  large  amounts  of  computer 
resources.  With  traditional  finite  element  analysis,  a 
structure  only  had  to  be  analyzed  once.  With  stochastic 
finite  element  analysis,  a  structure  might  have  to  be 
analyzed  hundreds,  even  thousands  of  times. 
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Problem  Statement 


Stochastic  Finite  Element  Analysis  (SFEA)  is  a  method  of 
finite  element  analysis  wherein  the  ultimate  product  of  the 
analysis  is  not  the  computed  stress  a  given  load  creates 
within  a  structure,  but  rather,  a  probability  that  the  struc¬ 
ture  will  successfully  withstand  a  given  load. 

Finite  Element  Analysis  (FEA)  is  a  tool.  Its  purpose  is  to 
simulate  real  world  conditions,  utilizing  a  mathematical 
model,  bolstered  up  by  underlying  assumptions  about  that 
model.  For  example,  we  can  generate  a  finite  element  mesh 
representing  a  wing.  The  model  consists  of  various  types  of 
quad  elements,  shear  panels  and  connecting  rods  to  model  the 
skin  surfaces,  ribs  and  spars  of  the  simulated  wing.  Then, 
for  purposes  of  the  analysis,  we  assume  the  wings  material 
properties  are  constant  and  that  the  loading  to  which  it  is 
subjected  represents  reality. 

These  models  may  be  very  good,  perhaps  the  result  of  years  of 
data  collection  and  study.  However,  when  the  model  and  its 
assumptions  are  coded  into  an  input  deck  for  the  finite 
element  analysis,  only  one  set  of  operational  conditions  is 
represented,  not  the  continuum  of  changing  conditions  present 
in  reality.  Moreover,  because  input  values  may  be  averages 
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or  modified  by  factors  of  safety,  the  simulation  of  reality 
may  be  inaccurate. 

Nevertheless,  the  results  of  the  finite  element  method  are 
generally  acceptable.  Buildings  don't  fall  down  often  and 
airplanes  operate  for  long  periods  without  their  wings 
tearing  off.  However,  when  reality  has  not  been  modeled 
closely  enough,  or  when  the  answers  are  not  precise  enough, 
the  need  to  depart  from  traditional  FEA  emerges.  Traditional 
finite  element  analysis  allows  neither  the  model  not  its 
assumptions  to  change.  SFEA  allows  the  designer  to  capture 
the  previously  discussed  modeling  and  structural  uncertain¬ 
ties  in  a  manner  more  closely  approximating  reality.  For 
example,  the  designer  knows  the  wind  loading  on  a  wing  is  not 
static  but  varies  by  small  gusts  and  buffets.  Perhaps  tests 
have  disclosed  that  a  certain  percentage  of  connections  will 
be  flawed  or  there  is  a  chance  certain  structural  elements 
may  be  weakened  or  destroyed  by  external  forces  during  opera¬ 
tion.  These  can  now  be  accounted  for. 

In  SFEA,  the  design  engineer  allows  these  parameters  to  vary 
according  to  known  -'or  assumed)  probability  distribution 
functions.  However,  because  the  internal  workings  of  finite 
element  programs  require  the  input  to  be  composed  of  deter¬ 
ministic  functions,  we  are  led  to  the  fundamental  difference 
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between  implementing  stochastic  and  traditional  finite 
element  models. 


Traditionally,  a  finite  element  analysis  is  run  only  once, 
the  output  of  stresses  and  strains  are  checked,  and  if 
needed,  adjustments  are  made  to  the  model  or  its  assumptions. 
The  process  may  be  repeated  to  check  the  adjustments,  but 
typically,  only  a  few  iterations  are  performed.  In  SFEA,  the 
process  is  repeated  a  large  number  of  times,  utilizing  Monte 
Carlo  simulation  (3:274)  to  model  the  uncertainties  in 
material  and  structural  properties.  Within  each  simulation 
run,  the  analysis  is  deterministic;  that  is,  it  is  performed 
exactly  like  any  other  finite  element  analysis.  Once  each 
run  is  completed,  the  results  are  stored  away  in  a  cumulative 
data  file  until  enough  data  has  been  collected  to  allow 
statistically  valid  calculations.  The  output  of  these  stan¬ 
dard  statistical  calculations  is  a  probability  that  an 
element  or  assembly  of  elements  will  not  exceed  some  pre¬ 


viously  defined  criterion. 


Methods  and  Results 


The  finite  element  model  used  in  this  demonstration  of  the 
method  of  SFEA  was  a  model  of  the  Intermediate  Complexity 
Wing,  (see  Figure  1)  designed  and  provided  by  the  Structures 
Division,  Wright  Research  and  Development  Center,  Wright- 
Patterson  Air  Force  Base,  OH.  It  depicts  the  wing  of  a  small, 
modern,  high  performance  fighter  aircraft,  somewhat  resembling 
an  F-5.  The  loading  is  aerodynamic  and  simulates  supersonic 
flight.  The  load  was  scaled  to  bring  the  most  highly  stressed 
elements  right  to  the  edge  of  maximum  allowable  stress.  This 
effectively  removed  any  excessive  design  by  reducing  the  fac¬ 
tor  of  safety  to  1.0.  This  hypothetical  wing's  purpose  is  to 
provide  engineers  with  a  reliable,  benchmark  model  to  help  in 
formulating  techniques  and  theories  and  for  use  in  compar¬ 
isons.  Its  intermediate  complexity  refers  to  the  size  of  the 
elemental  mesh.  It  is  fine  enough  to  yield  good  analysis 
results,  but  is  coarse  enough  to  be  reasonably  economical  in 
terms  of  computer  time. 

In  this  demonstration  of  technique,  the  analysis  was  oriented 
toward  evaluating  the  survivability  of  the  wing,  assuming  it 
has  suffered  the  loss  of  one  structura.1  element.  This  repre¬ 
sents  the  wing  suffering  damage  from  a  non-specific  external 
source.  This  could  mean  a  bird  strike,  being  hit  by  a  large 
caliber  projectile,  such  as  a  bullet,  or  any  other  type  of 
damage  which  would  destroy  the  structural  integrity  of  the 
element  during  flight.  Although  the  model  consists  of  158 
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elements,  only  119  were  considered  susceptible  to  damag  ;55 
shear  panels,  (Figure  2),  2  triangular  elements  and  62  quad 
elements,  (Figure  3)).  This  is  because  39  of  the  elements  are 
rod  elements  (Figure  4)  whose  function  is  to  hold  the  upper 
and  lower  wing  surfaces  apart.  In  reality,  ribs  and  spars 
would  perform  this  function,  but  in  the  model,  these  ribs  and 
spars  are  represented  by  shear  panels  incapable  of  resisting 
axial  stress.  Thus,  the  model  contains  members  not  present  in 
reality,  but  needed  to  make  the  model  computationally  correct. 

For  this  demonstration,  two  quantities  were  selected  as  random 
variables.  The  first  random  variable  was  the  structural  ele¬ 
ment  to  be  damaged.  Damage  was  simulated  by  a  reduction  in 
the  modulus  of  elasticity  of  the  randomly  selected  element  by 
a  factor  of  100;  it  could  no  longer  carry  any  significant 
stress.  The  other  random  variable  was  the  aerodynamic  loading 
at  each  node.  This  load  was  assumed  to  vary  according  to  a 
Gaussian,  or  normal,  distribution. 

At  the  beginning  of  each  analysis,  a  FORTRAN  coded  prepro¬ 
cessor,  called  PRE.FOR,  was  run  to  assemble,  from  an 
unaltered,  original  data  deck,  the  input  data  deck  needed  for 
the  Monte  Carlo  simulation.  The  preprocessor  generated  a 
random  number  corresponding  to  an  element  to  be  damaged.  As 
an  additional  evaluation  tool,  the  wing  was  divided  into 
approximate  thirds:  root,  mid,  and  tip  (see  Figures  1,  2,  3, 
and  4  for  divisions.)  The  wing  taken  as  a  whole,  that  is, 
damage  not  restricted  to  any  particular  area,  is  referred  to 
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as  the  combined  wing.  The  purpose  of  dividing  the  wing  into 
sections  was  to  allow  insight  into  possible  interrelations 
between  damage  in  one  element  and  stresses  in  another,  perhaps 
some  distance  away.  As  it  turned  out,  (and  is  discussed 
later) ,  there  indeed  was  a  significant  interelement  stress 
relation.  For  example,  damage  in  the  tip  section  resulted  in 
stress  failures  in  even  the  root  elements. 

For  the  first  one  third  of  the  runs,  damage  was  allowed  to 
occur  in  the  root  section  only.  If  the  randomly  selected 
element  was  not  within  the  root  section,  it  was  ignored  and 
another  was  picked.  Because  all  elements  in  a  finite  element 
model  are  uniquely  numbered,  this  discrimination  process  was 
simple.  Essentially,  a  uniform  random  number  generator 
(8:195-199)  produced  a  random  integer  between  1  and  119  (cor¬ 
responding  to  the  number  of  damageable  elements.)  This  number 
was  compared  to  a  list  of  elements  contained  within  the  root 
section.  Likewise,  for  the  second  series  of  runs,  only  the 
elements  within  the  mid  section  were  allowed  to  be  damaged. 

The  final  third  of  the  runs  had  damage  restricted  to  the  tip 
section.  (See  Appendix  D  for  more  detail  about  the  selection 
process . ) 

The  likelihood  that  any  one  element  type  was  selected  for 
damage  was  based  on  the  surface  that  particular  element  type 
presented  normal  to  the  surface  of  the  wing.  The  surface  area 
of  the  wing  in  the  "z”  or  "up”  direction  was  about  1925  square 
units.  The  area  of  the  front  of  the  wing  (the  leading  edge 
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shear  panels)  was  about  180  square  units.  The  surface  area  of 
the  shear  elements  in  the  "z”  direction  was  about  5  square 
units.  Adding  these  areas  together,  then  taking  the  approx¬ 
imate  fraction  of  the  shear  panel  area  to  the  total  area 
yielded  the  approximate  ratio  of  one  to  ten.  Thus,  large, 
flat  panels,  the  quad  elements,  were  about  ten  times  more 
likely  to  be  "hit"  than  the  thin  ribs  and  spars,  (i.e.,  the 
shear  elements) .  Additionally,  the  quad  and  triangular  ele¬ 
ments  were  always  damaged  in  pairs — the  top  panel  and  the 
panel  immediately  below  it — simulating  the  way  a  projectile 
might  pass  through  both.  This  partitioning  of  damage  neglects 
several  damage  modes  deemed  unnecessarily  complex  for  this 
demonstration.  For  example,  a  projectile  might  penetrate  the 
wing  at  an  oblique  angle,  damaging  a  lower  quad  panel,  one  or 
more  shear  panels,  and  possibly  one  or  more  additional  quad 
panels  on  its  way  out.  Or,  a  bird  strike  to  the  leading  edge 
could  easily  damage  several  leading  edge  shear  and  quad 
elements.  Nevertheless,  algorithms  could  be  written  to 
accomodate  these  effects.  Frequency  and  effects  of  bird- 
strikes  are  well  studied  and  the  effects  of  oblique  angle 
projectile  damage  could  be  coded  to  the  level  of  detail 
required. 

Once  a  random  element  within  the  proper  section  was  selected, 
the  preprocessor  read  through  the  original  data  deck  looking 
for  the  selected  element's  data  cards.  Once  found,  the  pre¬ 
processor  generated  new  cards  giving  the  affected  element,  or 
elements,  reduced  load  carrying  capacity  by  reducing  the 
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modulus  of  elasticity.  Next,  the  nodal  force  cards  were  read 
in.  The  original  forces  for  each  node  were  considered  the 
mean  for  that  node  and  the  preprocessor  varied  the  forces 
about  that  mean  with  an  associated  coefficient  of  variation  of 
0.1  (9).  Finally,  the  preprocessor  completed  the  new  data 
deck  by  copying  the  rest  of  the  original  data  deck,  such  as 
grid  coordinates  and  constraints,  to  the  end  of  the  new  deck. 

Next,  this  input  data  deck  was  fed  into  the  finite  element 
analysis  program  called  ASTROS  (Automated  STRuctural  Optimi¬ 
zation  System.)  Again,  this  was  supplied  by  the  Structures 
Division  of  WRDC.  ASTROS  was  instructed  to  perform  stress  and 
displacement  calculations  for  all  elements  within  the  model. 
Ultimately,  the  displacement  data  was  discarded  as  not  being 
useful  within  the  scope  of  this  demonstration.  Because  of  a 
lack  of  meaningful  displacement-oriented  failure  criteria 
(other  than  that  which  could  be  obtained  through  stress-strain 
calculations) ,  any  failure  mode  based  on  excessive  displace¬ 
ment  would  have  been  arbitrary  and  meaningless. 

At  the  end  of  the  analysis,  another  FORTRAN  program,  called 
POSTl.FOR  (for  postprocessor,  level  1)  searched  the  ASTROS 
output  data  file.  It  looked  for  elemental  stresses,  combining 
them  according  to  the  Von  Mises-Hencky  Stress  Criterion 
(10:85) 

stress  =  s  =  (x^  +  (1) 

where : 

X  =  major  principal  stress 

y  =  minor  principal  stress 
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Maximum  allowable  stress  was  computed  this  same  way  using  55 
ksi  allowable  tensile  stress  and  45  ksi  allowable  compressive 
stress,  for  a  maximum  allowable  stress  of  50.745  ksi  (6) . 

These  elemental  stresses  were  written  to  the  end  of  a  cumula¬ 
tive  raw  elemental  stress  data  file  for  further  analysis. 

The  entire  process  was  further  automated  by  invoking  a  simple 
batch  file  to  run  sequentially  PRE.FOR,  ASTROS,  and  POSTl.FOR 
a  specified  number  of  times,  with  an  average  run  time  for  the 
sequence  of  three  programs  of  about  65  cpu  seconds  (150 
seconds,  real  time.)  These  programs  were  run  on  a  Digital 
Equipment  Corporation  VAXstation  II  GPX  minicomputer, 
operating  at  3  mips. 

For  baseline  calculations,  the  loop  was  repeated  two  thousand 
times  for  each  of  the  three  wing  sections.  Then,  a  second 
postprocessor,  called  P0ST2.F0R  was  invoked  to  change  the  file 
of  raw  elemental  stress  data  into  statistical  data.  Computed 
for  each  element  was  its  mean  stress,  the  standard  deviation 
of  the  stress,  and  the  probability  that  it  would  be  stressed 
beyond  the  maximum  allowable  stress,  given  random  damage  some¬ 
where  on  the  wing.  The  stress  distributions  were  assumed  to 
be  normally  distributed  random  variables  thus  allowing  the 
probability  of  failure  to  be  approximated  by  (1:932) 

p^  =  1  -  .5(1  +  c^x  +  c^x^  +  c^x'^  )  (2) 

where:  x  =  (V  ”><)/(r 

•t*  =  maximum  allowable  stress 

»<=  mean  elemental  stress 

O'  =  standard  deviation  of  stress 

c  =  0.196854  c  =  0.115194 

c;  =  0.000344  ci  =  0.019527 

3  4 
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Once  these  elemental  probabilities  were  obtained,  they  were 
plotted  out  with  SURFER  (4),  a  three-dimensional  graphics 
program.  The  resulting  graphs  show  the  estimated  probability 
of  structural  failure  at  an  elemental  level  given  that  random 
damage  has  occurred  somewhere  in  the  wing  (Figure  5) .  Figures 
6,  7,  and  8  show  computed  probability  of  structural  failure 
given  that  damage  has  occurred  somewhere  in  the  root,  mid  and 
tip  sections  respectively.  Note  that  in  this  model,  damage 
even  to  the  tip  section  results  in  increased  probability  of 
failure  in  the  root  section.  Figure  8.  This  is  because  damage 
to  any  element  results  in  the  inability  of  that  element  to 
carry  its  share  of  nodal  loading.  Since  that  load  must  be 
picked  up  by  the  surrounding  elements,  its  effect  is  trans¬ 
mitted  through  moment  resisting  joints  quite  a  distance  away. 
Because  some  of  the  root  elements  were  highly  stressed  anyway, 
the  additional  stress  transmitted  to  them  was  sufficient  to 
push  them  past  the  failure  criterion. 

The  next  step  was  to  investigate  the  use  of  SFEA  as  a  tool  to 
increase  the  reliability,  or  in  this  case,  the  survivability 
of  a  structure.  Baseline  failure  probabilities  were  assembled 
using  2000  runs  per  section  of  the  wing,  for  a  total  of  6000 
runs.  Initially,  elemental  failure  probabilities  were  not 
known,  of  course,  but  were  estimated  to  be  in  the  range  of 
0.25  to  0.50,  given  the  highly  stressed  initial  conditions  in 
many  elements,  even  without  damage. 


The  following  equation  for  estimating  error  in  Monte  Carlo 
simulations  was  used  to  select  the  number  of  runs  to  establish 
the  baseline  (3:291-292) 

%  error  =  200  [(1  -  )/ (np^) ] (3) 

where: 

=  elemental  failure  probability 
=  number  of  iterations 

The  results  of  this  estimator  are  plotted  out  in  Figure  9  for 
several  selected  failure  probabilities.  Since  an  error 
between  5%  and  10%  was  desired,  this  led  to  a  range  of  number 
of  runs  from  400  to  4800.  Two  thousand  was  selected  as  a  good 
blend  of  accuracy  and  economy. 

Once  baseline  statistics  were  assembled,  the  two  quad  elements 
with  the  lowest  probability  of  surviving  random  damage  (ele¬ 
ments  63  and  64)  were  strengthened  by  increasing  their  thick¬ 
nesses  by  1/64  inch.  Then  the  entire  analysis  was  run  again. 
However,  this  time  through,  each  section  was  analyzed  only  200 
times,  instead  of  2000.  This  decreased  number  of  runs  meant 
higher  error,  up  to  a  maximum  of  100%  for  a  failure  proba¬ 
bility  of  0.02,  but  that  increase  would  be  offset  by  a  corres¬ 
ponding  decrease  in  the  failure  probabilities  through  the 
design  modifications.  Thus  even  a  100%  error  with  a  failure 
probability  of  0.02  still  yields  a  small  failure  probability 
range  (0-0.04) . 
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Percent  Error 


The  results  of  this  series  of  runs  were  analyzed.  The 
unexpected  result  was  that  while  the  additional  thickness  in 
the  two  strengthened  elements  raised  the  likelihood  they  would 
survive  a  hit,  the  additional  stiffness  imparted  by  the  extra 
thickness  caused  new  failures  in  near  by  elements.  Figure  10 
graphs  the  failure  probabilities  of  the  combined  wing  (all 
three  sections  susceptible  to  damage)  for  the  unimproved  wing 
and  the  wing  with  elements  63  and  64  strengthened.  Note  the 
very  low  failure  probability  in  elements  63  and  64  (right 
corner  of  graph)  and  the  slight  increase  in  failure  proba¬ 
bilities  in  elements  45,  46,  53,  54,  55,  and  56  (the  peak  next 
to  63  &  64) .  This  increase  in  failure  probabilities  is  more 
evident  when  the  wing  is  sectioned  off.  Figure  11  shows  a 
comparison  of  the  original  vs.  the  strengthened  wing  for 
damage  confined  to  the  root  area.  Not  much  change  is  evident, 
but  in  Figures  12  and  13,  showing  comparisons  of  the  mid  and 
tip  sections,  the  increase  in  failure  probabilities  is  ob¬ 
vious.  Effects  of  stress  telegraphing  are  plain,  especially 
with  tip  damage.  The  newly  strengthened  elements  are  failure- 
free  while  the  adjacent  element  has  been  pushed  into  failure. 
This  result  means  that  strengthening  parts  of  a  structure  can 
actually  lower  its  resistance  to  failure  under  certain  condi¬ 
tions,  such  when  damage  is  confined  to  a  certain  area.  The 
next  iteration  led  to  strengthening  these  neighboring  ele¬ 
ments,  another  round  of  600  total  runs,  followed  by  the  proba¬ 
bility  calculations.  This  time,  about  twenty  neighboring 
elements  were  strengthened.  A  trend  was  emerging.  While 
adding  strength  to  selected  elements  still  resulted  in 
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neighboring  elements  experiencing  more  failures,  the  trend  was 
flattening  out,  resulting  in  more  elements  being  more  sur- 
vivable.  Figure  14  compares  the  initial  improvement  (round 
one)  combined  wing  with  the  secondary  improvement  (round  two) 
combined  wing.  Overall,  the  failure  probabilities  have  de¬ 
creased.  Figure  15  compares  the  root  section  of  the  round  one 
and  round  two  wing.  Again  there  is  a  general  decrease  in 
failure  probabilities.  Figure  16  compares  the  mid  section  of 
the  two  wings.  While  failure  probabilities  are  still  rela¬ 
tively  high  in  the  mid  section  of  the  round  two  wing,  wing¬ 
breaking  stress  is  no  longer  being  shifted  into  other  sec¬ 
tions.  Figure  17  compares  the  tip  sections  of  the  two  wings. 
Note  that  even  though  failure  probabilities  increased  for  one 
element  in  the  mid  section  of  the  round  two  wing,  there  were 
no  longer  any  root  section  failures. 

By  the  end  of  the  third  round  of  strengthening,  most  elements 
had  reliabilities  of  over  ninety  percent,  and  those  that 
didn't  were  only  down  to  about  85%.  Note,  in  Figure  18  comp¬ 
aring  the  combined  wing  of  round  two  with  that  of  round  three, 
the  general  improvement  in  failure  probabilities.  Improvement 
is  evident  when  comparing  each  of  the  three  wing  subsections 
with  the  round  two  version,  in  Figures  19,  20,  and  21. 
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FIGURE  10.  Combined  Wing  Failure  Probabilities;  Original  vs 
Round  One  Improvements 
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Pr  o  tr  o  b  I  f  i  t  y 


The  final  step  was  to  trim  some  weight  from  those  elements 
whose  survivability  was  never  in  question.  Obviously,  these 
elements  were  in  a  lightly  stressed  state  as  their  relia¬ 
bilities  were  always  1.0.  With  this  final  adjustment,  the 
survivability  of  almost  every  element  was  greater  than  90%,  as 
opposed  to  the  original  design,  wherein  many  elements  had 
reliabilities  in  the  60-70%  range.  Interestingly,  trimming 
the  weight  from  these  lightly  stressed  elements  helped  improve 
the  reliability  of  more  heavily  stresses  neighboring  elements, 
in  much  the  same  way  that  adding  weight  to  heavily  stressed 
elements  in  round  one  decreased  the  reliability  of  adjacent 
elements.  On  the  downside,  however,  this  weight  savings  did 
result  in  some  of  the  tip  elements'  failure  probabilities 
increasing.  For  purposes  of  this  demonstration,  this  final 
step  results  in  a  wing  referred  to  as  "best”.  This  is  not 
meant  to  imply  that  this  is  as  reliable  as  the  wing  can  be, 
but  only  that  it  was  as  good  as  this  demonstration  made  it. 

Figure  22  shows  the  element  numbers  of  those  elements  whose 
strength  was  affected  in  each  round  of  design  modification. 

All  incremental  changes  were  1/64  inch  in  thickness.  Numbers 
in  paraentheses  represent  decreases  in  thickness  to  save 
weight.  The  cost  of  this  increase  in  survivability  was  esti¬ 
mated  at  a  mere  4%  increase  in  weight.  (See  Appendix  E  for 
more  details  of  this  estimation.)  This  was  estimated  by 
computing  the  solid  volume  of  the  structural  material  of  the 
original  wing  model  and  the  new  structural  solid  volume 
imparted  by  the  design  modifications. 
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LIST  OF  STRENGTHENED  ELEMENTS 


Round  One 


Round 


45, 

46 

19, 

20 

(  5, 

6) 

53, 

54 

21, 

22 

(  7, 

8) 

55, 

56 

27, 

28 

(  9, 

10) 

29, 

30 

(11, 

12) 

31, 

32 

(13, 

14) 

33, 

34 

(15, 

16) 

35, 

36 

(17, 

18) 

37, 

38 

(19, 

20) 

39, 

40 

(25, 

26) 

43, 

44 

(41, 

42) 

45, 

46 

51, 

52 

53, 

54 

57, 

58 

Parentheses  indicate  elements  lightened 


FIGURE  22.  Elements  Improved,  Round  by  Round 


In  contrast,  to  get  about  the  same  increase  in  survivability 
using  the  factor  of  safety  approach,  a  safety  factor  of  1.25 
had  to  be  applied  to  get  reliabilities  in  approximately  the 
same  range  as  the  best  wing  obtained  by  the  stochastic  method. 
It  may  be  interesting  to  note  the  usual  range  of  safety 
factor  used  for  structural  aluminum  is  1.65-2.34,  with  an 
average  about  2.0  (7:510).  While  the  safety  factor  for  a 
highly  weight-controlled  structure,  such  as  an  aircraft,  could 
be  expected  to  have  a  fairly  low  safety  factor,  it  is  not 
normally  below  1.5;  unmanned  structures,  such  as  missiles  can 
go  as  low  as  1.25  (6).  Figure  23  shows  the  combined  wing 
comparison  between  the  best  wing  obtained  stochastically  and 
the  wing  improved  by  the  1.25  factor  of  safety.  Overall,  the 
stochastic  wing  has  better  reliabilities.  Figure  24  and  25 
show  the  root  and  mid  section  comparisons.  The  stochastically 
best  wing  has  low  failure  probabilities  in  the  damaged  root 
and  mid  sections  along  with  a  few  very  low  failure  probabil¬ 
ities  in  the  tip  area,  as  shown  by  the  wavy  lines  at  the  tip. 
The  factor  of  safety  wing  had  higher  failure  probabilities  in 
the  damaged  sections,  but  essentially  zero  failure  proba¬ 
bilities  outside  the  damaged  sections.  Only  when  damage  was 
confined  to  the  tip  section.  Figure  26,  do  we  see  the  reverse 
of  this.  The  factor  of  safety  wing  had  very  low  failure 
probabilities  all  across  the  wing,  whereas  the  stochastically 
best  wing  had  tip  failures  rates  as  high  as  0.1.  Note,  how¬ 
ever,  a  safety  factor  of  1.25  means  that  weight  was  increased 
by  25% !  The  stochastic  method  of  analysis  resulted  in  a  more 
reliable  structure,  with  a  weight  savings  of  more  than  20% ! 
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The  final  step  in  this  demonstration  involved  comparing  two 
different  optimization  schemes.  One  scheme  used  the  opti¬ 
mizing  capability  of  ASTROS.  In  an  iterative  procedure, 

ASTROS  adjusts  the  thickness  of  the  model's  structural  members 
to  exactly  resist  the  applied  loads  (which  are  adjusted  to  the 
desired  factor  of  safety.)  ASTROS  is  one  of  the  very  few 
finite  element  programs  available  with  this  built-in  function; 
virtually  all  other  programs  rely  on  the  designers  ability  to 
interpret  the  output  data  and  update  the  model  accordingly. 
Thus,  ASTROS  represents  the  state-of-the-art  in  optimization 
codes . 

The  other  method  was  to  continue  the  stochastic  process 
already  started,  using  a  a  short  program,  P0ST3.F0R,  to  update 
the  elements'  thicknesses  based  on  failure  probabilities. 

First,  ASTROS  was  instructed  to  optimize  the  model  to  a  safety 
factor  of  1.25.  The  estimated  weight  gain  over  the  original 
model  was  about  22%,  while  the  failure  probabilities  compared 
to  the  stochastically  best  wing  (at  4%  weight  gain)  are  shown 
in  Figure  27.  Next,  ASTROS  optimized  to  a  safety  factor  of 
1.5,  with  a  resultant  weight  gain  estimated  at  about  30%  over 
the  original  wing.  The  failure  probabilities  of  this  wing 
were  essentially  zero  out  to  the  fourth  and  fifth  decimal 
place.  Then  the  previous  stochastic  programs  were  run  under  a 
new  batch  file  along  with  the  POST3.FOR  program  acting  as  a 
crude  optimizer.  (The  term  crude  is  applied  as  POST3.FOR  did 
not  allow  for  any  lightening  of  the  elements — it  would  only 
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strengthen  them.)  Essentially,  the  stochastic  optimizer 
worked  like  this.  At  the  end  of  the  standard  600  runs, 
P0ST2.F0R  was  called  to  generate  elemental  failure  proba¬ 
bilities.  Then,  P0ST3.F0R  would  scan  that  probability  list 
looking  for  "large”  failure  rates.  (What  constituted  large 
changed  as  the  failure  probabilities  converged  to  zero.  For 
the  first  couple  of  runs,  large  was  0.1;  later,  it  was  0.001. 
The  number  could  be  changed  within  P0ST3.F0R  while  the 
optimizer  was  running.)  It  would  then  increment  the  thickness 
of  the  failure  prone  elements,  and  the  whole  procedure  would 
start  again,  for  ten  iterations.  At  this  point,  the 
stochastically  optimized  wing's  failure  probabilities  very 
closely  matched  the  ASTROS  optimized  wing  with  the  safety 
factor  of  1.5.  The  weight  gain  for  this  stochastically  opti¬ 
mized  wing  was  about  25%.  Thus,  for  essentially  the  same 
reliability,  the  stochastic  method  of  design  utilizing  a  crude 
optimizing  scheme  gave  a  design  10%  lighter  than  that  produced 
by  one  of  the  best  optimization  programs  available. 
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Conclusions 


Stochastic  Finite  Element  Analysis  is  a  tool  which  may  be  used 
by  designers  to  attempt  to  improve  the  reliability  or  the 
weight  of  a  structure.  Its  use  involves  some  limitations  and 
drawbacks  in  comparison  to  traditional  methods  of  design  using 
finite  element  analysis  and  factors  of  safety.  Principally, 
the  limitations  are  in  the  designer's  ability  to  adequately 
capture  probability  distribution  functions  of  potential  fail¬ 
ure  modes  and  then  to  encode  them  into  a  computer  digestable 
algorithm.  The  principal  drawback  is  the  potential  for  using 
very  large  amounts  of  computer  time. 

In  this  demonstration,  the  limitations  were  overcome  by  using 
simple  probability  distribution  functions  to  model  the 
hazards.  Random  variables  were  restricted  to  just  two:  damage 
location  and  nodal  loading.  The  random  variable  of  damage  was 
restricted  to  only  one  or  two  structural  elements  that  were 
easy  to  select.  More  complex  modes  of  damage  would  involve 
commensurately  more  complex  programming,  especially  if  mul¬ 
tiple  damage  sites  and  cascading  damage  effects  are  consi¬ 
dered.  The  random  variable  associated  with  nodal  loading  was 
confined  to  an  easily  modeled  normal  distribution.  Only  one 
loading  condition  was  considered  here,  that  of  supersonic 
flight.  Many  other  high  stress  loading  conditions  exist,  such 
as  high-g  maneuvers,  heavily  loaded  munitions  pylons,  even 
shock  wave  loading  from  nearby  explosions.  Conceivably,  a 
wing  optimized  for  one  loading  condition  could  be  a  terrible 
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performer  for  another.  However,  since  structural  design  is  a 
study  in  compromise,  a  solution,  representing  improvement  over 
the  traditional  approach,  should  be  almost  always  achievable. 

The  drawback  of  using  vast  amounts  of  computer  time  could  be 
insurmountable.  Baseline  analysis  of  the  Intermediate  Complexity  Wing 
took  a  week  of  CPU  time  on  a  three  mip  minicomputer,  almost  16 
days  real  time,  and  several  weeks  human  time.  Manipulation  of 
the  resulting  enormous  data  files  was  very  slow.  Subsequent 
600  iteration  runs  took  eighteen  hours,  real  time. This  simple 
demonstration  used  at  least  fifteen  of  them,  not  including 
false  starts,  mistakes,  and  debugging  runs.  A  very  large, 
complex  structure,  such  as  a  bomber  wing,  utilizing  thousands 
of  nodes,  could  take  a  week  for  just  one  run.  Access,  there¬ 
fore,  to  extremely  powerful,  inexpensive  computers  is  essen¬ 
tial  to  successfully  and  economically  apply  the  method  of 
stoch?.stic  finite  element  analysis  to  any  large,  complex 
structure. 

Nevertheless,  given  reasonably  accurate  hazard  probabilities 
and  access  to  sufficient  computer  resources,  this  method  can 
be  a  valuable  to  in  providing  a  lighter,  more  reliable 
structure.  Figures  27  and  28  trace  the  reliability  growth  of 
a  wing  originally  designed  using  traditional  finite  element 
analysis  methods.  As  elements  of  the  wing  were  selectively 
strengthened,  the  wing's  failure  probabilities  were  system¬ 
atically  reduced  to  a  fairly  uniform  low  rate.  Figure  29 
compares  the  stochastically  strengthened  wing  to  the  wing 
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strengthened  through  the  factor  of  safety.  While  both  are 
highly  reliable,  the  factor  of  safety  wing  is  more  than  20% 
heavier  then  the  stochastic  wing. 

When  ASTROS,  a  state-of-the-art  finite  element  program,  was 
called  upon  to  produce  an  optimized  wing  based  on  a  safety 
factor  of  1.5,  it  yielded  a  wing  of  extremely  high  relia¬ 
bility.  The  stochastic  algorithms  developed  for  this  demon¬ 
stration  also  yielded  a  wing  of  similar  reliability,  but  with 
a  weight  savings  of  10%  over  the  ASTROS  optimized  wing. 
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Recommendations 


Several  recomendations  were  derived  from  this  study: 

(1)  The  entire  scope  could  be  expanded  in  a  follow  on  study. 
The  model  doesn't  have  to  be  any  larger,  but  additional  ele¬ 
ment  types,  or  more  complex  geometry,  might  provide  additional 
insight  to  the  utility  of  the  method. 

(2)  Introduce  additional  levels  of  complexity  into  the  random 
variables.  For  example,  have  more  than  one  loading  condition 
available,  and  allow  the  way  elements  are  damaged  to  expand. 

(3)  Introduce  additional  random  variables  into  the  analysis. 
Let  a  variable  representing  fatigue  cracking  enter  the  pic¬ 
ture,  for  example. 

(4)  Obtain  a  finite  element  code  whose  output  could  be  more 
controlled  than  ASTROS.  Its  output  always  contained  a  lot  of 
extraneous  material,  such  as  page  headers,  footers,  and  the 
like  that  had  to  be  sifted  through  to  get  to  the  real  data. 
That  added  to  the  programming  complexity,  as  well  as  added 
time  to  each  run. 

(5)  Try  to  find  a  higher  speed  computer.  The  computer  used 
in  this  study  was  adequate  for  this  work,  but  would  not  be 
fast  enough  in  a  significantly  more  complex  project. 
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Program  Listings 
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REAL  LOC,MEAN 
CHARACTER  REST* 6 4 
COMMON  INIT,  I SEED 
INIT— 1 
C 

C  Note  the  input  file  is  broken  up  into  small  files  for  easier  handl 
C 

OPEN  ( 6, FILE- ' velmnt. dat STATUS- 'OLD' ) 

OPEN  (7,FILE-'icw.inp' ,STATUS-'OLD' ) 

OPEN  (8,FILE-'iseed' ,  STATUS- ' OLD ' ) 

OPEN  (9,FILE-'vforce.dat' ,STATUS-'OLD' ) 

OPEN  ( 10, FILE-' control. dat' ,STATUS-'OLD' ) 

OPEN  ( ll,FILE-'vgrid.dat' ,STATUS-'OLD' ) 

OPEN  (12,FILE-'vprop.dat' ,STATUS-'OLD' ) 

OPEN  ( 14, FILE-' broke. elem' , STATUS- 'OLD' ) 

C 

READ(8,1000)ISEED 

C 

5  READ( 10, 500, END-10)  REST 

WRITE(7,500)  rest 
GO  TO  5 
C 

10  CALL  FORCE 

C 

C  *URAND*  returns  uniform  r.n.  to  localize  damage  in  wing, 

C  based  approximately  on  percentages  of  areas  normal  to  wing 
C 

15  LOC-URAND  (I SEED) 

LOC-LOC*100 

C 

IF  (LOC.LT.90.0)  CALL  QUAD(ISEED) 

IF  (90.0. LE.LOC.AND.LOC.lt. 99.0)  CALL  SHEAR(ISEED) 

IF  (99.0.LE.LOC.AND.ISEED.GT.4000)  then 
CALL  TRIANG(ISEED) 

elseif( 99.0. le.loc.and.iseed.lt. 4 001)  then 
go  to  15 
END  IF 
C 
C 

1000  FORMAT(I7) 

500  FORMAT(A64) 

C 

C 

STOP 

END 

C 

C 

SUBROUTINE  TRIANG( ISEED) 

CHARACTER  TYPE*8,  POINTS*32,  REST*64 
INTEGER  ELEM,  MAT,  PROP 


52 


Feb  15  07:28  1989  pre.for  Page  2 


C 

WRITE(14, 2000)1, 2 

C  Looks  for  triangular  element  and  sub  in  new  mat  property 
20  READ  (6, 1500, END-35)  TYPE, ELEM, MAT, POINTS 

IF  (LGE(TYPE, 'CTRMEM' ) )  THEN 
MAT- 2 
END  IF 

WRITE  (7,1500)  TYPE, ELEM, MAT, POINTS 
GO  TO  20 
C 
C 

35  READ( 12, 500, END-40)  REST 
WRITE(7,500)  REST 
GO  TO  35 
C 

40  READ (11, 500, END-45) REST 

WRITE(7,500)REST 
GO  TO  40 
C 
C 

2000  F0RMAT(2I5) 

1500  F0RMAT(A8,2I8,A32) 

1000  F0RMAT(I3,5X,I4) 

500  FORMAT(A64) 

C 

C 

45  RETURN 
END 
C 
C 

SUBROUTINE  SHEAR(ISEED) 

CHARACTER  TYPE*8,  POINTS*32,  REST*64 
INTEGER  ELEN,  MAT,  SHRPNL,  MATER,  COUNTER 
TEMP-URAND( I SEED) 

SHRPNL-TEMP*55+65 

C  Identifies  one  of  55  shear  panels  and  adds  65  to  it  to 

C 

C 

IF( ISEED.lt. 2001)CALL  ROOT(SHRPNL) 

IF( ISEED.GT.2000  .AND.  ISEED. LT. 4001 )CALL  MID(SHRPNL) 
IF(ISEED.GT.4001)CALL  TIP( SHRPNL) 

C 

C  Looks  for  shear  panels  in  general  and  specifically  for 
C  shear  panel  matching  location  given  by  SHRPNL 
C 

WRITE( 14 , 2000 )SHRPNL 
C 

40  READ  (6, 1500, end-50)  TYPE, ELEM, MAT, POINTS 

IF(LGE(TYPE, 'C8HEAR' ) .AND. LLE( TYPE , 'CSHEAR' ) .AND. 
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+  ELEM.EQ.SHRPNL)  THEN 

MATER-MAT 
MAT- 4 4 
END  IF 

WRITE  (7,1500)  TYPE, ELEM, MAT, POINTS 
GO  TO  40 
C 

C  Looks  for  PSHEAR  card  associated  with  damaged  panel 
50  READ  (12, 1500, END-55)  TYPE, ELEM, MAT, POINTS 

IF(LGE(TYPE, 'PSHEAR' ) . AND. LLE ( TYPE , 'PSHEAR' ) .AND. 

+  ELEM. EQ. MATER)  THEN 

WRITE (7,1500) TYPE ,44,3, POINTS 
END  IF 

WRITE  (7,1500)  TYPE, ELEM, MAT, POINTS 
GO  TO  50 
C 

55  READ( 11, 500, END-60 )REST 

WRITE(7,500)REST 
GO  TO  55 
C 

2000  F0RMAT(I5) 

1500  FORMAT  (A8,2l8,A32) 

1000  FORMAT  (I3,5X,I4) 

500  FORMAT  (A64) 

C 

C 

60  RETURN 
END 
C 
C 

SUBROUTINE  QUAD(ISEED) 

CHARACTER  TYPE*8,  POINTS*32,  REST*64 
INTEGER  ELEM , MAT , MEMBl , MEMB2 , NUMB 
C  Sets  up  the  two  specific  quadmems  to  be  reset 

57  TNUMB-URAND(ISEED) 

NUMB-TNUMB*30 

MEMB2-(NUMB-i-2)*2 

MEMB1-MEMB2-1 

IF(ISEED.LE.2000  .AND.  MEMBl. GE. 41)  THEN 
GO  TO  58 

ELSEIF( ISEED.GE.2001  .AND.  ISEED.LE. 4000  .AND.  MEMBl. GE. 25  .AND. 
•«-MEMBl.LE.40)THEN 
GO  TO  58 

ELSEIF(ISEED.GE.4001  .AND.  MEMBl. LE. 23)  THEN 
GO  TO  58 
END  IF 
GO  TO  57 

C  Looks  for  specific  cquad  member  and  the  one  below  it. 

C 

58  WRITE(14,2000)MEMB1,NEMB2,  ISEED 
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60  READ  (6, 1500, END-70)  TYPE, ELEM, MAT, POINTS 

IF(LLE(TyPE, 'CQDMEMl  ') -AND. ELEM . EQ . MEMBl . OR. ELEM . EQ .MEMB2 )  THEN 
C 

MATER-MAT 
MAT- 4 3 
C 

END  IF 

WRITE  (7,1500)  TYPE, ELEM, MAT, POINTS 
GO  TO  60 

C  Reset  their  property  cards 
C 

70  READ  (12, 1500, END-75)  TYPE, ELEM, MAT, POINTS 

IF(LLE(TYPE, 'PQDMEMl  ') -AND. ELEM. EQ .MATER )  THEN 
WRITE  (7,1500)  TYPE, 43, 3, POINTS 
END  IF 

WRITE  (7,1500)  TYPE, ELEM, MAT, POINTS 
GO  TO  70 
C 

75  READ( 11, 500, END-80)  REST 
WRITE(7,500)  REST 
GO  TO  75 
C 

2000  F0RMAT(2I5,I5) 

1500  F0RMAT(A8,2I8,A32) 

1000  F0RMAT(I3,5X,I4) 

500  FORMAT (A6 4) 

C 

80  RETURN 
END 
C 
C 

SUBROUTINE  FORCE 
C 

CHARACTER  FRCE*32,  VECTOR* 2 4 
REAL  MAG,  MEAN 
C 

80  READ  (9, 1500, END-85)  FRCE, MAG, VECTOR 

SDEV-0.1*MAG 
Y-GAUSS(  MAG,  SDEV  ) 

MAG-Y 

WRITE(7,1500)  FRCE, MAG, VECTOR 
GO  TO  80 

1500  FORMAT(A32,F8.2,A24) 

C 

C 

85  RETURN 

END 
C 
C 

SUBROUTINE  ROOT( SHRPNL ) 
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INTEGER  SHRPNL 

90  IF( SHRPNL. GE. 89  .AND.  SHRPNL . LE . 96 ) RETURN 
IF(SHRPNL.GE.101  .AND.  SHRPNL . LE . 103 ) RETURN 
IF(SHRPNL.GE.109  .AND.  SHRPNL . LE . Ill ) RETURN 
IF(WHRPNL.GE.117  .AND.  SHRPNL. LE . 119 ) RETURN 
TEMP-URAND( I SEED) 

SHRPNL-TEMP*  5  5+6  5 
GO  TO  90 
END 
C 

SUBROUTINE  MID (SHRPNL) 

INTEGER  SHRPNL 

91  IF( SHRPNL. GE. 78  .AND.  SHRPNL . LE . 88 ) RETURN 
IF(SHRPNA.GE.99  .AND.  SHRPNL . LE . 100 ) RETURN 
IF(SHRPNL.GE.107  .AND.  SHRPNL . LE . 108 ) RETURN 
IF(SHRPNL.GE.114  .AND.  SHRPNL . LE . 115 ) RETURN 
TEMP-URAND ( I SEED ) 

SHRPNL-TEMP* 55+6 5 

GO  TO  91 

END 

SUBROUTINE  TIP (SHRPNL) 

INTEGER  SHRPNL 

92  IF(SHRPNL.GE.65  .AND.  SHRPNL. LE. 76 ) RETURN 
IF(SHRPNL.GE.97  .AND.  SHRPNL . LE . 98 ) RETURN 
IF(SHRPNL.GE.105  .and.  SHRPNL. LE. 106) return 
IF(SHRPNL.GE.113  .AND.  SHRPNL . LE . 114 ) RETURN 
TEMP-URAND ( I SEED ) 

SHRPNL-TEMP*55+64 
GO  TO  92 
END 
C 

REAL  FUNCTION  GAUSS  (  MEAN,  STD  ) 

— GAUSS— ———— ——————————————— 

REAL  MEAN 

COMMON  INIT 


DATA  I SET  /O/ 

♦INITIALIZATION  OF  -URAND-  FOR  THE  RANDOM  NUMBER  SEQUENCE. 
IF  (  ISET  .EQ.  0  )  THEN 
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#THE  FOLLOWING  COMPUTATIONS  PRODUCE  TWO  RANDOM  DEVIATES.  RETURN 
ONE  AND  SAVE  THE  OTHER  FOR  THE  NEXT  CALL.  USE  ODDCAL  AS  A  CUE 
FOR  WHAT  TO  DO  NEXT  TIME  IN. 

#PICK  TWO  UNIFORM  DEVIATES  IN  THE  SQUARE  EXTENDING 
FROM  -1  TO  +1. 

100  VI  -  2.*URAND(ISEED)  -  1. 

V2  -  2.*URAND{ ISEED)  -  1. 

#SEE  IF  VI  AND  V2  ARE  IN  THE  UNIT  CIRCLE.  IF  NOT,  TRY  AGAIN. 
R  -  Vl**2  +  V2**2 
IF  (  R  .GE.  1.  )  GO  TO  100 

#MAKE  THE  BOX-MULLER  TRANSFORMATION  TO  PRODUCE  TWO  NORMAL 
DEVIATES  WITH  MEAN  ZERO  AND  STANDARD  DEVIATION  ONE. 

TEMP  -  SQRT(-2.*LOG(R)/R) 

ENORM  -  VI* TEMP 

ONORM  «  V2*TEMP 

GAUSS  -  MEAN  +  STD* ONORM 

ISET  -  1 

ELSE 

#USE  ENORM  FROM  PREVIOUS  CALL  TO  COMPUTE  GAUSS. 

GAUSS  •  MEAN  +  STD*ENORM 
ISET  -  0 
END  IF 

RETURN 

END 

REAL  FUNCTION  URAND  (ISEED) 

URAND—— « 
COMMON  I NIT 


INTEGER  Al 

/ 

A2 

9 

A3 

INTEGER  Cl 

9 

C2 

9 

C3 

INTEGER  Ml 

9 

H2 

9 

M3 

PARAMETER 

(Ml  - 

259200,  Al 

m 

7141, 

Cl 

-  54773) 

PARAMETER 

(M2  - 

134456,  A2 

m 

8121, 

C2 

-  28411) 

PARAMETER 

(M3  - 

243000,  A3 

m 

4561, 

C3 

-  51349) 

PARAMETER 

(NSHUF 

-  30) 

REAL  UDEV(NSHUF) 

RMl  -  l./Ml 
RM2  -  1./M2 
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C  #INITIALIZE  ON  FIRST  CALL. 

IF  (  INIT  .LT.  0  )  THEN 

OPEN(8,FILE-'iseed' , STATUS- 'OLD' ) 

REWIND(8) 

READ(8,115)  ISEED 
REWIND! 8) 

ISEED-IO+ISEED 
WRITE(8,115)  ISEED 
115  FORMAT! 17) 

II  -  MOD!ABS!Cl-ISEED) ,  Ml) 

11  -  M0D!A1*I1+C1,  Ml) 

12  -  MOD! II,  M2) 

II  -  MOD(Al*Il+Cl,  Ml) 

13  -  MOD!Il,  M3) 

tINITIALIZE  THE  SHUFFLING  VECTOR  FOR  RANDOM  SEQUENCE. 

DO  11  I  -  1,  NSHUF 

11  -  MOD!Al*Il+Cl,  Ml) 

12  -  M0D!A2*I2+C2,  M2) 

UDEV!I)  -  ! REAL! II)  +  REAL ! 1 2 ) *RM2 ) *RMl 
11  CONTINUE 

INIT  -  1 
ENDIF 

tGENERATE  THE  NEXT  NUMBER  IN  EACH  INTEGER  SEQUENCE  II,  12,  I3. 

11  -  MOD!Al*Il+Cl,  Ml) 

12  -  M0D!A2*I2+C2,  M2) 

13  -  M0D!A3*I3+C3,  M3) 

#USE  THE  THIRD  INTEGER  SEQUENCE  TO  GET  A  NUMBER  BETWEEN 
1  AND  NSHUF. 

I  -  1  +  !NSHUF*I3)  /  M3 

#RETURN  THAT  VECTOR  ENTRY. 

URAND  -  UDEV!I) 

#REFILL  SHUFFLING  VECTOR  WHERE  SAMPLE  HAS  BEEN  USED. 

UDEV!I)  -  ! REAL! II)  +  REAL! I2 ) *RM2 ) *RMl 

RETURN 

END 
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CHARACTER  INPUT*80,  DATA* 48,  TARGET*! 0,  BUSTED* 10 
INTEGER  J,  N,  LOOP,  POINT 

REAL  STRESS(158) , DISP ( 88 } , PRESSl , PRESS2 , MAXSTR 
REAL  MAXD1,MAXD3,MAXD5,MAXD7,MAXD9 
DATA  J/l/ 

DATA  N/0/ 

DATA  NSET/0/ 

DATA  LOOP/0/ 

DATA  POINT/0/ 

DATA  LOOPD/0/ 

PARAMETER  ( MAXSTR-55750 . 0 ) 

PARAMETER  (MAXDl-12.0,  MAXD3-12.75,  MAXD5-13.5) 
PARAMETER  ( MAXD7-14 . 25 ,  MAXD9-15.0) 

PARAMETER  (NNODE-88) 

OPEN(10,FILE-'icw,out' ,  STATUS- 'OLD' ) 
OPEN(ll,FILE-'Stress.dat' ,STATUS-'OLD' ) 
OPEN(12,FILE-'displ.dat' , STATUS- ' OLD ' ) 
OPEN(14,FILE-' iseed' ,STATUS-'OLD' ) 

OPEN (15, FILE- 'damage' , STATUS- 'OLD' ) 
OPEN(20,FILE-'broke.elem' ,STATUS-'OLD' ) 

C 

READ  (14, 1000) ISEED 
READ  (20, 1200) BUSTED 
C 
C 

100  READ  (10, 2000) INPUT 
C 

C  Looks  for  node  heading  on  displacement  page 
POINT-INDEX( INPUT, 'POINT  ID.') 

200  IF(POINT.GT.0)THEN 
GO  TO  225 
ELSE 

GO  TO  100 
END  IF 
C 

225  LOOPD-LOOPD+1 

GO  TO  (300,300)LOOPD 
C 

250  READ  (10,2000)INPUT 

C  Looks  for  element  heading  on  stress  pages 
POINT-INDEX( INPUT, 'ELEMENT' ) 

IF( POINT. GT.O)  THEN 
GO  TO  275 

Er  ^ 

GO  TO  250 
ENDIF 

275  LOQP-LOOP+1 
C 

C  Recycles  to  look  for  various  elemental  stress  pages 

GO  TO  (400,600,800,900)  LOOP 
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C 

C  Computes  displacement  magnitude  for  each  element  and  writes  to  file 
C 

300  READ(10,3000,ERR-100)  J , DISPLl , DISPL2 , DISPL3 

DISP( J)-SQRT(DISPL1**2+DISPL2**2+DISPL3**2) 

IF( J.EQ.NNODE)THEN 
GO  TO  250 
ELSE 

GO  TO  300 
END  IF 
C 

C  This  loop  reads  and  processes  the  quad  member  stresses 
C  Following  acts  as  line  feed  to  get  to  lines  of  data 
400  READ(10,2000)  INPUT 

C  Following  ERR  statement  jumps  to  second  page  of  quad  stresses 
500  READ(10,5000,ERR-250)  M, PRESSl , PRESS2 

C  Following  computes  Mises-Hinkley  Stress  Criterion,  writes  it  to  fill 
STRESS ( M ) -SQRT ( PRESSl **2+PRESS2  *  *2-PRESSl *PRESS2 ) 

IF(M.EQ.64 )THEN 
GO  TO  250 
else 

GO  TO  500 
endi  f 
C 

C  Following  acts  as  a  line  feed 
C 

600  READ( 10, 2000) INPUT 
C  Following  loop  reads  rod  stresses 
625  READ( 10,6000 )  M,  STRESS(M),  STRESS(M+1) 

C 

IF  (M,LT.156)  THEN 
GO  TO  625 
C 

ELSE 

READ(10,8000)  M,STRESS(M) 

GO  TO  250 
C 

ENDIF 

C 

800  READ{10,2000)  INPUT 

C 

C  Following  reads  shear  panel  stresses 
C 

825  READ(10,7000)  M, STRESS (M) , STRESS (M+1 ) 

IF(M.LT.117)  THEN 
GO  TO  825 
ELSE 

READ(10,7500)  M,STRESS(M) 

GO  TO  250 
C 
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END  IF 
C 

GO  TO  825 
C 

C  Following  reads  triangular  panel  stresses 
C 

900  READ(10,2000)  INPUT 

C 

925  READ(10,5000,err-10)  M,PRESS1 , PRESS2 

C 

STRESS(M)-SQRT(PRESS1**2+PRESS2**2-PRESS1*PRESS2) 
lF(M.eq.2)go  to  10 
C 

GO  TO  925 
C 

C  Following  puts  pointers  at  end  of  cumulative  data  files 
C  and  writes  new  data  lines. 

C 

10  READ( 12, 2000) INPUT 

C 

c  GO  TO  10 

C 

cO  WRITE(12,4000)  (DISP(K),  K-1,88) 

C 

C 

30  READ( 11, 2000, END-40)  INPUT 

C 

GO  TO  30 
C 

40  WRITE(11,9000)  (STRESS(I),  1-1,158) 

C 

50  READ (15, 2000, END-55) INPUT 

GO  TO  50 
C 

55  DO  60  1-1,158 

I F  (  STRESS ( I ) . GT . MAXSTR ) THEN 
NSET  -  1 
END  IF 

60  CONTINUE 

C 

IF(NSET.EQ.l )  then 

go  to  62 

else 

WRITE ( 15,2500 ) I  SEED, BUSTED 
END  IF 
C 

62  DO  65  1-1,158 

I F ( STRESS ( I ) . GT . MAXSTR ) THEN 

WRITE( 15,1500 ) ISEED, BUSTED, I , STRESS (  I ) 

END  IF 
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65 

c 

c 

c 

c 

c 

c 

c 

c 

c 

cO 

c 

1000 

1200 

1500 

1600 

2000 

2500 

3000 

4000 

5000 

6000 

7000 

7500 

8000 

9000 

C 


CONTINUE 

i£( dispd  )  .gt  .maxdl )  write (15, 1600 )  iseed,  1  ,disp(  1 ) 
i£(disp(  3  )  . gt  .maxd3  )  writed5,1600  )  iseed,  3  ,disp(  3  ) 
i£(disp(  5)  .gt.inaxd5)  write(  15, 1600) iseed,  5,disp(  5) 
i£(disp(7)  gt.maxd7)  write(15,1600)iseed,7,disp(7) 
i£(disp( 9 ) . gt .maxd9 )  write (15, 1600 ) iseed, 9,disp( 9 ) 

DO  70  1-1,9, 2 

IF(DISP(I) .GT.MAXDIS(I) )THEN 
WRITE(15,1600)ISEED,1,DISP(I) 

ENDIF 

CONTINUE 

FORMAT (17) 

FORMAT (AlO ) 

FORMAT( 'RUN' ,15, '  BRKN  ELEMS ' , AlO , '  FLD  ELEM' , 1 5 , '  STR-' 
FORMAT('RUN  #  ',15,'  NODE  #  ',15,'  DISPLACEMENT-  ',E14. 
FORMAT (A80 ) 

FORMAT( 'Run' , 15,  '  Broken  Elements:  ',A10,'  No  Failures' 
F0RMAT(11X,I7,15X,3E14.6) 

F0RMAT(88E14.6) 

F0RMAT(I9,64X,2E14.6) 

F0RMAT(I13,E14.6,45X,E14.6) 

FORMAT( I14,E16 .6,45X,E16.6) 

F0RMAT(I14,E16.6) 

F0RMAT(I13,E14.6) 

F0RMAT(158E14.6) 

STOP 

END 


,E14. 

6) 

) 
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real  stress( 6000, 158) ,sdev(158) , mean (158) ,sum(158) 
real  sum2 ( 158 ) ,prob( 158 ) 

C  Loops  is  the  total  number  of  iterations  you've  done  w/  Astros 

C  Remember  to  update  this  number  to  reflect  the  total  number  of  runs 

C  to  be  evaluated 

parameter  (loops  -  200 , sigmax-52500 ) 

parameter  (cl-. 196854, c2-. 115194 , c 3-. 000344, c4-. 019527) 

open( 10, file- 'stress', status-' old' ) 
open( 20, file- 'output' , status- 'old' ) 
open( 30 , file- 'prob.dat ' , status-' old' ) 

do  100  i-1, loops 

read( 10,1000, enJ-100 ) ( stress( i , j ) , j-1 , 158 ) 

100  continue 

do  300  j-1, 158 
do  200  i-1, loops 

sum( j )-stress( i , j )+sum( j ) 
sum2 ( j )-stress(i,j ) **2+sum2( j ) 

200  continue 
300  continue 

do  400  j-1, 158 

mean( j )-sum( j )/loops 

sdev( j )-sqrt( sum2( j )/loops-mean( j )**2) 


400  continue 

do  500  j-1, 158 

write (20, 2000)j ,mean( j ) , sdev( j ) 

500  continue 

do  600  m-1,158 

X- ( s i gmax-mean ( m ) ) /sde v ( m ) 

prob(m)-l-. 5* ( i+ci*x+c2*x**2+c3*x**3+c4*x**4 ) **-4 
wri te( 30 , 3000 )m,prob(m) 

600  continue 


1000  format(158El4.6) 

2000  format ( 'mean  in  element' ,I4, '  is',El4.6,'  etd  dev  is',El4.6) 

3000  format(i5,f8.3) 

stop 

end 
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c  This  program  transforms  a  list  of  probabilties 

c  into  an  input  file,  coordinate  vs.  probabilty, 

c  for  input  into  a  3-D  graphics  program,  additional 

c  points  are  generated  to  give  the  graph  a  smoother  surface. 

character  misc*64 

real  x(88) ,y(88) ,xbar(1000) ,ybar(1000) ,prob(158) 

open( 10 , file- 'vgrid.dat' , status- ' old' ) 
open( 20, file- 'velmnt.dat' , status- 'old' ) 
open( 30, f ile-' prob.dat ' , status- 'old' ) 
open( 40 , file- 'wing. grid' , status- 'old' ) 
c  Gets  past  first  five  lines  to  get  to  the  element 
c  and  its  corresponding  node  numbers. 

do  5  m-1,5 
read( 10 , 5000 )misc 
5  continue 

c  reads  in  element  and  its  node  numbers 
do  10  m-1,88 
read(10,1000)n,x(n) ,y(n) 

10  continue 

do  15  m-1,119 

read(  30, 3000)prob(m) 

15  continue 

c  finds  center  of  triangular  element 

read ( 20 , 2000)n,i,j,k,l 
xbar(n)-(x( i )+x( j )+x(k) )/3 
ybar(n)-(y( i )+y( j )+y(k) )/3 

read( 20, 5000 )misc 

c  following  loop  locates  four  corners  of  quad  elements 
c  and  assigns  coordinates  to  them 

do  20  m-3,63,2 

read ( 20 ,2000)n,i,j,k,l 

ex-  abs( x( i )-x( j ) ) 
wy-  abs(y( i )-y(k) ) 

xbar ( m)-( x( i )+x( j )+x( k )+x{ 1 ) )/4 
xbar(m+100)-  x(i)+.l*ex 
xbar(m+200)-  x(j)-.l*ex 
xbar(m+300)-  x(k)-.l*ex 
xbar(m+400)-  x^)  +  .l*ex 

ybar (m)-(y( i )+y( j )+y(k)+y( 1 ) )/4 
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ybar(m+100;-  y(i)-.l*wy 
ybar(m+200)-  y(j)-.l*wy 
ybar(m+300)«  y(k)+.l*wy 
ybar(m+400)-  y(l)+.l*wy 

read( 20 , 5000  )  misc 

20  continue 

c  following  loop  locates  the  line  of  each  shear  element 

c  and  gives  coordinates  to  several  points.  figures  out 

c  which  direction  the  element  is  pointing  to  make  sure 

c  the  coordinates  follow  the  line  correctly, 

do  30  m-1,55 

read ( 20 ,2000)n,i,j,k,l 

ex-  abs( x( i )-x( k  )  ) 
wy-  abs(y(i)-y(k) ) 

xbar (m+500 )«x( i )  + . 01*x( i  ) 
xbar (m+600 )-x( k )  +  . 01*x( k ) 
xbar (m+700 )-( X ( i )+x( k ) )/2 

if  (x( i ) .gt.x(k)  )  then 

xbar(m+800)»x(k)+.25*ex 

xbar(m+900)-x(k)+.75*ex 

else 

xbar(m+800)-x( i )+.25*ex 
xbar(m+900)-x( i )+.75*ex 

endif 

ybar  (.’+500  )-y(  i  )  +  .  01*y(  i  ) 
yba  r( m+600) -y(k)  +  .01*y(k) 
ybar(m+700)-(y(i)+y(k) )/2 

if  (y(i) .gt.y(k) )  then 

ybar (m+800)-y( k )+.75*wy 
ybar (m+900 )-y( k )+ . 25*wy 

else 

yba  r(m+8C0)-y(i)  +  .25  *wy 
ybar (m+900 )-y( i )+. 75*wy 

endif 
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30  continue 

c  writes  each  coordinate  and  its  associated  prob  to  file 
write( 40, 4000 )xbar ( 1 ) ,ybar (1) ,prob( 1 ) 

do  400  m«3,63,2 

do  450  n-0,400,100 

write( 40 , 4000 )  xbar(m+n),  ybar(m+n),  prob(m) 

450  continue 

400  continue 

do  550  m«l,55 
do  500  n-500,900,100 

wri  te  (  40 , 4000  )xbar  (m+n)  ,ybar  (m+n)  ,prob(in+64  ) 

500  continue 
550  continue 


1000  format(8x,i8,8x,2f8.3) 

2000  format(8x,i8,8x,4i8) 
3000  format(5x,f8.3) 

4000  format(f8.3,f8.3,f8.3) 

5000  format(a64) 

stop 

end 
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c  This  program  provides  a  means  to  use  pre.for, 

c  ASTROS,  postl.for,  run.x,  and  post2.for  in  an 

c  optimizing  scheme,  with  optimize. x  as  the 

c  controller.  Note  this  program  does  not  allow 

c  an  element  to  be  reduced  in  weight. 

character  type*8 

open( 10, f ile-'prob.dat' ,status-'old' ) 
open ( 20 , file-' vprop.dat' , status- 'old' ) 
open( 30, f ile-' vp.dat ' , status- 'old' ) 
open (40, file- 'St ress.dat' , status- ' old' ) 
open( 50 , f ile- ' iseed' , status-' old' ) 

10  read( 20 , 1000 , end-100 ) type,n,mat, thick 

read( 10,2000 )num,prob 
x-abs( 1 . 0-prob) 

if (x.ge.0.005)thick-thick+0.005 
write( 30, 3000) type, n, mat, thick 
go  to  10 

100  write ( 40 , 4000 ) 

write(50, 5000)1 

1000  format(a8,2i8,f8.4) 

2000  format(i5,f8.3) 

3000  format(a8,2i8,f8.4) 

4000  format( '  '  ) 

5000  format(i5) 
end 


Feb  25  00:16  1989  run.x 


# 

set  n-0 

while  ($n  !-  400) 
pre.x 

astros  icw.inp 
postl  .X 
@  n++ 
end 


Feb  25  00:16  1989  optimize. x 


i 

set  m-0 

while  ($m  !-  5) 
run.x 

p6  0  0 .  X 
po6t3  .X 

cp  vp.dat  vprop.dat 
@  m++ 
end 
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ASSIGN  DATABASE  ICW  TEST  NEW  DELETE 
SOLUTION 

TITLE-  INTERMEDIATE  COMPLEXITY  WING 

ANALYZE 

BOUNDARY  SPC  -  1 
STATICS  (MECH  -  1) 

PRINT  DISP-ALL,  STRESS-ALL 


END 

BEGIN  BULK 
FORCE 

1 

3 

15.17 

0.0000 

0.0000 

FORCE 

1 

4 

15.26 

0.0000 

0.0000 

FORCE 

1 

5 

-1538.43 

1.0000 

0.0000 

FORCE 

1 

5 

-3924.44 

0.0000 

1.0000 

FORCE 

1 

5 

570.04 

0.0000 

0.0000 

FORCE 

1 

6 

1393.27 

1.0000 

0.0000 

FORCE 

1 

6 

4481.66 

0.0000 

1.0000 

FORCE 

1 

6 

618.14 

0.0000 

0.0000 

FORCE 

1 

7 

48.87 

0.0000 

0.0000 

FORCE 

1 

8 

50.47 

0.0000 

0.0000 

FORCE 

1 

9 

-5165.61 

1.0000 

0.0000 

FORCE 

1 

9 

-4908.70 

0.0000 

1.0000 

FORCE 

1 

9 

490.31 

0.0000 

0.0000 

FORCE 

1 

10 

4636.29 

1.0000 

0.0000 

FORCE 

1 

10 

5067.86 

0.0000 

1.0000 

FORCE 

1 

10 

566.79 

0.0000 

0.0000 

FORCE 

1 

1 

106.29 

1.0000 

0.0000 

FORCE 

1 

1 

-4164.07 

0.0000 

1.0000 

FORCE 

1 

1 

489.42 

0.0000 

0.0000 

FORCE 

1 

2 

-120.12 

1.0000 

0.0000 

FORCE 

1 

2 

3445.47 

0.0000 

1.0000 

FORCE 

1 

2 

424.33 

0.0000 

0.0000 

FORCE 

1 

11 

86.57 

0.0000 

0.0000 

FORCE 

1 

12 

101.77 

0.0000 

0.0000 

FORCE 

1 

13 

131.35 

0.0000 

0.0000 

FORCE 

1 

14 

133.78 

0.0000 

0.0000 

FORCE 

1 

15 

149.66 

0.0000 

0.0000 

FORCE 

1 

16 

119.26 

0.0000 

0.0000 

FORCE 

1 

17 

-3099.10 

1.0000 

0.0000 

FORCE 

1 

17 

1034.29 

0.0000 

1.0000 

FORCE 

1 

17 

532.08 

0.0000 

0.0000 

FORCE 

1 

18 

2693.56 

1.0000 

0.0000 

FORCE 

1 

18 

-1333.79 

0.0000 

1.0000 

FORCE 

1 

18 

510.58 

0.0000 

0.0000 

FORCE 

1 

19 

1224.43 

1.0000 

0.0000 

FORCE 

1 

19 

-457.83 

0.0000 

1.0000 

FORCE 

1 

19 

414.60 

0.0000 

0.0000 

FORCE 

1 

20 

-1367.42 

1.0000 

0.0000 

FORCE 

1 

20 

527.02 

0.0000 

1.0000 

FORCE 

1 

20 

343.27 

0.0000 

0.0000 

FORCE 

1 

21 

169.08 

0.0000 

0.0000 

1.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

1.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000 

1.0000 
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FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 


177.12 

0.0000 

0.0000 

1.0000 

172.57 

0.0000 

0.0000 

1.0000 

149.96 

0.0000 

0.0000 

1.0000 

181.42 

0.0000 

0.0000 

1.0000 

200.56 

0.0000 

0.0000 

1.0000 

-2252.94 

1.0000 

0.0000 

0.0000 

963.52 

0.0000 

1.0000 

0.0000 

501.60 

0.0000 

0.0000 

1.0000 

1665.66 

1.0000 

0.0000 

0.0000 

-864.64 

0.0000 

1.0000 

0.0000 

420.19 

0.0000 

0.0000 

1.0000 

833.74 

1.0000 

0.0000 

0.0000 

-330.92 

0.0000 

1.0000 

0.0000 

323.85 

0.0000 

0.0000 

1.0000 

-893.03 

1.0000 

0.0000 

0.0000 

394.44 

0.0000 

1.0000 

0.0000 

290.28 

0.0000 

0.0000 

1.0000 

181.44 

0.0000 

0,0000 

1.0000 

161.47 

0.0000 

0.0000 

1.0000 

185.69 

0.0000 

0.0000 

1.0000 

207.31 

0.0000 

0.0000 

1.0000 

198.26 

0.0000 

0.0000 

1.0000 

172.47 

0.0000 

0.0000 

1.0000 

-2347.27 

1.0000 

0.0000 

0.0000 

893.76 

0.0000 

1.0000 

0.0000 

561.65 

0.0000 

0.0000 

1.0000 

2261.60 

1.0000 

0.0000 

0.0000 

-941.98 

0.0000 

1.0000 

0.0000 

479.01 

0.0000 

0.0000 

1.0000 

945.21 

1.0000 

0.0000 

0.0000 

-421.50 

0.0000 

1.0000 

0.0000 

390.19 

0.0000 

0.0000 

1.0000 

-778.56 

1.0000 

0.0000 

0.0000 

353.86 

0.0000 

1.0000 

0.0000 

347.96 

0.0000 

0.0000 

1.0000 

206.87 

0.0000 

0.0000 

1.0000 

181.70 

0.0000 

0.0000 

1.0000 

189.96 

0.0000 

0.0000 

1.0000 

237.12 

0.0000 

0.0000 

1.0000 

200.79 

0.0000 

0.0000 

1.0000 

206.23 

0.0000 

0.0000 

1.0000 

2182.20 

1.0000 

0.0000 

0.0000 

968.95 

0.0000 

1.0000 

0.0000 

536.20 

0.0000 

0.0000 

1.0000 

2296.42 

1.0000 

0.0000 

0.0000 

-946.42 

0.0000 

1.0000 

0.0000 

494.70 

0.0000 

0.0000 

1.0000 

1103.36 

1.0000 

0.0000 

0.0000 

-452.79 

0.0000 

1.0000 

0.0000 

417,00 

0.0000 

0.0000 

1.0000 

FORCE 

1 

50 

-912.25 

1.0000 

0.0000 

0.0000 

FORCE 

1 

50 

462.44 

0.0000 

1.0000 

0.0000 

FORCE 

1 

50 

435.68 

0.0000 

0.0000 

1.0000 

FORCE 

1 

51 

234.88 

0.0000 

0.0000 

1.0000 

FORCE 

1 

52 

194.34 

0.0000 

0.0000 

1.0000 

FORCE 

1 

53 

167.27 

0.0000 

0.0000 

1.0000 

FORCE 

1 

54 

204.93 

0.0000 

0.0000 

1.0000 

FORCE 

1 

55 

223.26 

0.0000 

0.0000 

1.0000 

FORCE 

1 

56 

232.96 

0.0000 

0.0000 

1.0000 

FORCE 

1 

57 

-2246.54 

1.0000 

0.0000 

0.0000 

FORCE 

1 

57 

956.33 

0.0000 

1.0000 

0.0000 

FORCE 

1 

57 

612.16 

0.0000 

0.0000 

1.0000 

FORCE 

1 

58 

2202.02 

1.0000 

0.0000 

0.0000 

FORCE 

1 

58 

-989.04 

0.0000 

1.0000 

0.0000 

FORCE 

1 

58 

625.45 

0.0000 

0.0000 

1.0000 

FORCE 

1 

59 

1314.61 

1.0000 

0.0000 

0.0000 

FORCE 

1 

59 

-494.60 

0.0000 

1.0000 

0.0000 

FORCE 

1 

59 

515.95 

0.0000 

0.0000 

1.0000 

FORCE 

1 

60 

-1233.54 

1.0000 

0.0000 

0.0000 

FORCE 

1 

60 

502.27 

0.0000 

1.0000 

0.0000 

FORCE 

1 

60 

312.37 

0.0000 

0.0000 

1.0000 

FORCE 

1 

61 

230.44 

0.0000 

0.0000 

1.0000 

FORCE 

1 

62 

211.40 

0.0000 

0.0000 

1.0000 

FORCE 

1 

63 

212.10 

0.0000 

0.0000 

1.0000 

FORCE 

1 

64 

184.19 

0.0000 

0.0000 

1.0000 

FORCE 

1 

65 

208.58 

0.0000 

0.0000 

1.0000 

FORCE 

1 

66 

195.87 

0.0000 

0.0000 

1.0000 

FORCE 

1 

67 

-1676.00 

1.0000 

0.0000 

0.0000 

FORCE 

1 

67 

745.21 

0.0000 

1.0000 

0.0000 

FORCE 

1 

67 

383.15 

0.0000 

0.0000 

1.0000 

FORCE 

1 

68 

1418.60 

1.0000 

0.0000 

0.0000 

FORCE 

1 

68 

-507.59 

0.0000 

1.0000 

0.0000 

FORCE 

1 

68 

501.30 

0.0000 

0.0000 

1.0000 

FORCE 

1 

69 

1696.74 

1.0000 

0.0000 

0.0000 

FORCE 

1 

69 

-232.75 

0.0000 

1.0000 

0.0000 

FORCE 

1 

69 

553.35 

0.0000 

0.0000 

1.0000 

FORCE 

1 

70 

-1440.36 

1.0000 

0.0000 

0.0000 

FORCE 

1 

70 

246.59 

0.0000 

1.0000 

0.0000 

FORCE 

1 

70 

577.07 

0.0000 

0.0000 

1.0000 

FORCE 

1 

71 

240.01 

0.0000 

0.0000 

1.0000 

FORCE 

1 

72 

186.24 

0.0000 

0.0000 

1.0000 

FORCE 

1 

73 

208.65 

0.0000 

0.0000 

1.0000 

FORCE 

1 

74 

160.27 

0.0000 

0.0000 

1.0000 

FORCE 

1 

75 

180.50 

0.0000 

0.0000 

1.0000 

FORCE 

1 

76 

153.32 

0.0000 

0.0000 

1.0000 

FORCE 

1 

77 

-670.41 

1.0000 

0.0000 

0.0000 

FORCE 

1 

77 

145.01 

0.0000 

1.0000 

0.0000 

FORCE 

1 

77 

232.81 

0.0000 

0.0000 

1.0000 

FORCE 

1 

78 

795.53 

1.0000 

0.0000 

0.0000 

FORCE 

1 

78 

-133.64 

0.0000 

1.0000 

0.0000 
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FORCE 

1 

78 

240.06 

0.0000 

0.0000 

CTRMEM 

1 

1 

1 

3 

11 

CTRMEM 

2 

1 

2 

4 

12 

CQDMEMl 

3 

3 

3 

5 

13 

11 

CQDMEMl 

4 

3 

4 

6 

14 

12 

CQDMEMl 

5 

3 

5 

7 

15 

13 

CQDMEMl 

6 

3 

6 

8 

16 

14 

CQDMEMl 

7 

4 

7 

9 

17 

15 

CQDMEMl 

8 

4 

8 

10 

18 

16 

CQDMEMl 

9 

3 

1 

11 

21 

19 

CQDMEMl 

10 

3 

2 

12 

22 

20 

CQDMEMl 

11 

3 

11 

13 

23 

21 

CQDMEMl 

12 

3 

12 

14 

24 

22 

CQDMEMl 

13 

5 

13 

15 

25 

23 

CQDMEMl 

14 

5 

14 

16 

26 

24 

CQDMEMl 

15 

5 

15 

17 

27 

25 

CQDMEMl 

16 

5 

16 

18 

28 

26 

CQDMEMl 

17 

5 

19 

21 

31 

29 

CQDMEMl 

18 

5 

20 

22 

32 

30 

CQDMEMl 

19 

6 

21 

23 

33 

31 

CQDMEMl 

20 

6 

22 

24 

34 

32 

CQDMEMl 

21 

6 

23 

25 

35 

33 

CQDMEMl 

22 

6 

24 

26 

36 

34 

CQDMEMl 

23 

7 

25 

27 

37 

35 

CQDMEMl 

24 

7 

26 

28 

38 

36 

CQDMEMl 

25 

4 

29 

31 

41 

39 

CQDMEMl 

26 

4 

30 

32 

42 

40 

CQDMEMl 

27 

8 

31 

33 

43 

41 

CQDMEMl 

28 

8 

32 

34 

44 

42 

CQDMEMl 

29 

8 

33 

35 

45 

43 

CQDMEMl 

30 

8 

34 

36 

46 

44 

CQDMEMl 

31 

8 

35 

37 

47 

45 

CQDMEMl 

32 

8 

36 

38 

48 

46 

CQDMEMl 

33 

6 

39 

41 

51 

49 

CQDMEMl 

34 

6 

40 

42 

52 

50 

CQDMEMl 

35 

9 

41 

43 

53 

51 

CQDMEMl 

36 

9 

42 

44 

54 

52 

CQDMEMl 

37 

10 

43 

45 

55 

53 

CQDMEMl 

38 

10 

44 

46 

56 

54 

CQDMEMl 

39 

11 

45 

47 

57 

55 

CQDMEMl 

40 

11 

46 

48 

58 

56 

CQDMEMl 

41 

7 

49 

51 

61 

59 

CQDMEMl 

42 

7 

50 

52 

62 

60 

CQDMEMl 

43 

10 

51 

53 

63 

61 

CQDMEMl 

44 

10 

52 

54 

64 

62 

CQDMEMl 

45 

12 

53 

55 

65 

63 

CQDMEMl 

46 

12 

54 

56 

66 

64 

CQDMEMl 

47 

13 

55 

57 

67 

65 

CQDMEMl 

48 

13 

56 

58 

68 

66 

CQDMEMl 

49 

7 

59 

61 

71 

69 

1.0000 


72 


CQDMEMl 

50 

7 

60 

62 

72 

70 

CQDMEMl 

51 

11 

61 

63 

75 

71 

CQDMEMl 

52 

11 

62 

64 

74 

72 

CQDMEMl 

53 

14 

63 

65 

75 

73 

CQDMEMl 

54 

14 

64 

66 

76 

74 

CQDMEMl 

55 

15 

65 

67 

77 

75 

CQDMEMl 

56 

15 

66 

68 

78 

76 

CQDMEMl 

57 

6 

69 

71 

81 

79 

CQDMEMl 

58 

6 

70 

72 

82 

80 

CQDMEMl 

59 

16 

71 

73 

83 

81 

CQDMEMl 

60 

16 

72 

74 

84 

82 

CQDMEMl 

61 

17 

73 

75 

85 

83 

CQDMEMl 

6.. 

17 

74 

76 

86 

84 

CQDMEMl 

63 

18 

75 

77 

87 

85 

CQDMEMl 

64 

18 

76 

78 

88 

86 

CSHEAR 

65 

20 

1 

2 

4 

3 

CSHEAR 

66 

20 

3 

4 

6 

5 

CSHEAR 

67 

20 

5 

6 

8 

7 

CSHEAR 

68 

20 

7 

8 

10 

9 

CSHEAR 

69 

20 

1 

2 

12 

11 

CSHEAR 

70 

20 

11 

12 

14 

13 

CSHEAR 

71 

20 

13 

14 

16 

15 

CSHEAR 

72 

20 

15 

16 

18 

17 

CSHEAR 

73 

20 

19 

20 

22 

21 

CSHEAR 

74 

20 

21 

22 

24 

23 

CSHEAR 

75 

20 

23 

24 

26 

25 

CSHEAR 

76 

20 

25 

26 

28 

27 

CSHEAR 

77 

20 

29 

30 

32 

31 

CSHEAR 

78 

20 

31 

32 

34 

33 

CSHEAR 

79 

20 

33 

34 

36 

35 

CSHEAR 

80 

20 

35 

36 

38 

37 

CSHEAR 

81 

20 

39 

40 

42 

41 

CSHEAR 

82 

20 

41 

42 

44 

43 

CSHEAR 

83 

20 

43 

44 

46 

45 

CSHEAR 

84 

20 

45 

46 

48 

47 

CSHEAR 

85 

20 

49 

50 

52 

51 

CSHEAR 

86 

20 

51 

52 

54 

53 

CSHEAR 

87 

20 

53 

54 

56 

55 

CSHEAR 

88 

20 

55 

56 

58 

57 

CSHEAR 

89 

21 

59 

60 

62 

61 

CSHEAR 

90 

22 

61 

62 

64 

63 

CSHEAR 

91 

20 

63 

64 

66 

65 

CSHEAR 

92 

20 

65 

66 

68 

67 

CSHEAR 

93 

23 

69 

70 

72 

71 

CSHEAR 

94 

24 

71 

72 

74 

73 

CSHEAR 

95 

20 

73 

74 

76 

75 

CSHEAR 

96 

20 

75 

76 

78 

77 

CSHEAR 

97 

20 

1 

2 

20 

19 

CSHEAR 

98 

25 

19 

20 

30 

29 

CSHEAR 

99 

26 

29 

30 

40 

39 

73 


CSHEAR 

100 

27 

CSHEAR 

101 

28 

CSHEAR 

102 

29 

CSHEAR 

103 

30 

CSHEAR 

104 

31 

CSHEAR 

105 

32 

CSHEAR 

106 

33 

CSHEAR 

107 

34 

CSHEAR 

108 

35 

CSHEAR 

109 

36 

CSHEAR 

110 

37 

CSHEAR 

111 

38 

CSHEAR 

112 

20 

CSHEAR 

113 

44 

CSHEAR 

114 

39 

CSHEAR 

115 

40 

CSHEAR 

116 

35 

CSHEAR 

117 

36 

CSHEAR 

118 

41 

CSHEAR 

119 

42 

CROD 

120 

19 

CROD 

121 

19 

CROD 

122 

19 

CROD 

123 

19 

CROD 

124 

19 

CROD 

125 

19 

CROD 

126 

19 

CROD 

127 

19 

CROD 

128 

19 

CROD 

129 

19 

CROD 

130 

19 

CROD 

131 

19 

CROD 

132 

19 

CROD 

133 

19 

CROD 

134 

19 

CROD 

135 

19 

CROD 

136 

19 

CROD 

137 

19 

CROD 

138 

19 

CROD 

139 

19 

CROD 

140 

19 

CROD 

141 

19 

CROD 

142 

19 

CROD 

143 

19 

CROD 

144 

19 

CROD 

145 

19 

CROD 

146 

19 

CROD 

147 

19 

CROD 

148 

19 

CROD 

149 

19 

39 

40 

50 

49 

49 

50 

60 

59 

59 

60 

70 

69 

69 

70 

80 

79 

5 

6 

14 

13 

13 

14 

24 

23 

23 

24 

34 

33 

33 

34 

44 

43 

43 

44 

54 

53 

53 

54 

64 

63 

63 

64 

74 

73 

73 

74 

84 

83 

9 

10 

18 

17 

17 

18 

28 

27 

27 

28 

38 

37 

37 

38 

48 

47 

47 

48 

58 

57 

57 

58 

68 

67 

67 

68 

78 

77 

77 

78 

88 

87 

i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

74 


CROD 

150 

19 

61 

62 

CROD 

151 

19 

63 

64 

CROD 

152 

19 

65 

66 

CROD 

153 

19 

67 

68 

CROD 

154 

19 

69 

70 

CROD 

155 

19 

71 

72 

CROD 

156 

19 

73 

74 

CROD 

157, 

19 

75 

76 

CROD 

158 

19 

77 

78 

PTRMEM 

1 

1 

.0312 

PTRMEM 

2 

3 

.0312 

PQDMEMl 

3 

2 

.0312 

PQDMEMl 

4 

2 

.0416 

PQDMEMl 

5 

2 

.0364 

PQDMEMl 

6 

2 

.0624 

PQDMEMl 

7 

2 

.0520 

PQDMEMl 

8 

2 

.0884 

PQDMEMl 

9 

2 

.1040 

PQDMEMl 

10 

2 

.1040 

PQDMEMl 

11 

2 

.1144 

PQDMEMl 

12 

2 

.1875 

PQDMEMl 

13 

2 

.1040 

PQDMEMl 

14 

2 

.1144 

PQDMEMl 

15 

2 

.1563 

PQDMEMl 

16 

2 

.0884 

PQDMEMl 

17 

2 

.1248 

PQDMEMl 

18 

2 

.2031 

PROD 

19 

1 

.0200 

PSHEAR 

20 

1 

.0190 

PSHEAP. 

21 

1 

.0220 

PSHEAR 

22 

1 

.0210 

PSHEAR 

23 

1 

.0250 

PSHEAR 

24 

1 

.0240 

PSHEAR 

44 

3 

.0380 

PSHEAR 

25 

1 

.0380 

PSHEAR 

26 

1 

.0420 

PSHEAR 

27 

1 

.0480 

PSHEAR 

28 

1 

.0470 

PSHEAR 

29 

1 

.0390 

PSHEAR 

30 

1 

.0310 

PSHEAR 

31 

1 

.0280 

PSHEAR 

32 

1 

.0370 

PSHEAR 

33 

1 

.0500 

PSHEAR 

34 

1 

.0580 

PSHEAR 

35 

1 

.0650 

PSHEAR 

36 

1 

.0790 

PSHEAR 

37 

1 

.1010 

PSHEAR 

38 

1 

.1260 

PSHEAR 

39 

1 

.0440 

PSHEAR 

40 

1 

.0530 

75 


PSHEAR  41  1  .0920 

PSHEAR  42  1  .1000 

GRDSET  456 


MATl 

1  10.5E+06 

.300 

.000259 

MATl 

2  10.5E+06 

.300 

.000259 

MATl 

3  10.5E+04 

.300 

.000259 

SPCl 

1 

123 

79 

THRU  1 

B8 

GRID 

1 

63.500 

90.000 

1.125 

GRID 

2 

63.500 

90.000 

-1.125 

GRID 

3 

70.833 

90.000 

1.313 

GRID 

4 

70.833 

90.000 

-1.313 

GRID 

5 

78.167 

90.000 

1.500 

GRID 

6 

78.167 

90.000 

-1.500 

GRID 

7 

85.500 

90.000 

1.313 

GRID 

8 

85.500 

90.000 

-1.313 

GRID 

9 

92.833 

90.000 

1.125 

GRID 

10 

92.833 

90.000 

-1.125 

GRID 

11 

69.686 

87.471 

1.349 

GRID 

12 

69.686 

87.471 

-1.349 

GRID 

13 

76.097 

84.851 

1.586 

GRID 

14 

76.097 

84.851 

-1.586 

GRID 

15 

82.746 

82.133 

1.427 

GRID 

16 

82.746 

82.133 

-1.427 

GRID 

17 

89.647 

79.312 

1.259 

GRID 

18 

89.647 

79.312 

-1.259 

GRID 

19 

57.266 

77.669 

1.279 

GRID 

20 

57.266 

77.669 

-1.279 

GRID 

21 

63.992 

74.920 

1.532 

GRID 

22 

63.992 

74.920 

-1.532 

GRID 

23 

70.962 

72.071 

1.799 

GRID 

24 

70.962 

72.071 

-1.799 

GRID 

25 

78.191 

69.116 

1.617 

GRID 

26 

78.191 

69.116 

-1.617 

GRID 

27 

85.692 

66.050 

1.424 

GRID 

28 

85.692 

66.050 

-1.424 

GRID 

29 

51.032 

65.339 

1,433 

GRID 

30 

51.032 

65.339 

-1.433 

GRID 

31 

58.297 

62.369 

1.715 

GRID 

32 

58.297 

62.369 

-1.715 

GRID 

33 

65.826 

59.291 

2.012 

GRID 

34 

65.826 

59.291 

-2.012 

GRID 

35 

73.635 

56.100 

1.807 

GRID 

36 

73.635 

56.100 

-1.807 

GRID 

37 

81.738 

52.787 

1.590 

GRID 

38 

81.738 

52.787 

-1.590 

GRID 

39 

44.799 

53.008 

1.587 

GRID 

40 

44.799 

53,008 

-1.587 

GRID 

41 

52.603 

49.818 

1.898 

GRID 

42 

52.603 

49.818 

-1.898 

GRID 

43 

60.691 

46.512 

2.225 
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GRID 

44 

60.691  46.512 

GRID 

45 

69.079  43.083 

GRID 

46 

69.079  43.083 

GRID 

47 

77.784  39.525 

GRID 

48 

77.784  39.525 

GRID 

49 

38.565  40.678 

GRID 

50 

38.565  40.678 

GRID 

51 

46.908  37.267 

GRID 

52 

46.908  37.267 

GRID 

53 

55.555  33.732 

GRID 

54 

55.555  33.732 

GRID 

55 

64.523  30.067 

GRID 

56 

64.523  30.067 

GRID 

57 

73.830  26.262 

GRID 

58 

73.830  26.262 

GRID 

59 

32.331  28.347 

GRID 

60 

32.331  28.347 

GRID 

61 

41.214  24.716 

GRID 

62 

41.214  24.716 

GRID 

63 

50.420  20.953 

GRID 

64 

50.420  20.953 

GRID 

65 

59.967  17.050 

GRID 

66 

59.967  17.050 

GRID 

67 

69.876  13.000 

GRID 

68 

69.876  13.000 

GRID 

69 

25.166  14.173 

GRID 

70 

25.166  14.173 

GRID 

71 

35.583  12.304 

GRID 

72 

35.583  12.304 

GRID 

73 

46.181  10.403 

GRID 

74 

46.181  10.403 

GRID 

75 

56.964  8.469 

GRID 

76 

56.964  8.469 

GRID 

77 

67.938  6.500 

GRID 

78 

67.938  6.500 

GRID 

79 

18.000  0.000 

GRID 

80 

18.000  0.000 

GRID 

81 

30.000  0.000 

GRID 

82 

30.000  0.000 

GRID 

83 

42.000  0.000 

GRID 

84 

42.000  0.000 

GRID 

85 

54.000  0.000 

GRID 

86 

54.000  0.000 

GRID 

87 

66.000  0.000 

GRID 

88 

66.000  0.000 

ENDDATA 

-2.225 

1.997 

-1.997 

1.756 

-1.756 

1.742 

-1.742 

2.082 

-2.082 

2.438 

-2.438 

2.187 

-2.187 

1.922 

-1.922 

1.896 

-1.896 

2.265 

-2.265 

2.651 

-2.651 

2.376 

-2.376 

2.088 

-2.088 

2.073 

•2.073 

2.446 

■2.446 

2.827 

■2.827 

2.502 

•2.502 

2.169 

2.169 

2.250 

2.250 

2.625 

2.625 

3.000 

3.000 

2.625 

2.625 

2.250 

2.250 
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Sample  ASTROS  Output 
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AUTOMATED  STRUCTURAL  OPTIMIZATION  SYSTEM 


INITIAL  PRODUCTION  RELEASE  -  VERSION  2 

DEC  VAX/VMS  SERIES 
JUNE  1,  1988 


ASSIGN  DATABASE  ICW  RON  NEW  DELETE 


INTERMEDIATE  COMPLEXITY  WING 


8  0  U 

N  0  A  R  Y  CO 

N  0  1  T  I  0  N 

SUMMARY 

•  •  •  • 

statics/normal 

MOOES 

•  •  • 

••  STATICS 

MASS  MODES 

SAERO 

•••  flutter 

•  YES 

NO  NO 

NO 

•  NO 

ASTROS  VERSION  2 

8/29/88 

FOR  BOUNDARY 

C  0 

N  D  I 

T  I  0  N  1 

DYNAMIC  RESPONSE 

•  ••• 

MODAL 

DIRECT 

TRANS  FREQ  BLAST 

•  •• 

TRANS 

FREQ 

BLAST  •• 

NO  NO  NO 

• 

NO 

NO 

NO  • 

79 


^‘t.RMLUlAl£  COMriExlTY  WING 


DISPLACE 


POINT  ID. 

TYPE 

T1 

T2 

1 

C 

-1 .532566-01 

-3.384676-01 

2 

G 

1 .532566-01 

3.384676-01 

3 

G 

-1 .669216-01 

-3.643676-01 

4 

G 

1 .669216-01 

3.643876-01 

6 

0 

-2.115646-01 

-4.694166-01 

6 

G 

2.115646-01 

4.694166-01 

7 

G 

-1 .838646-01 

-3.952996-01 

6 

G 

1  .838646-01 

3.962996-01 

9 

G 

-1 .630576-01 

-3.999446-01 

10 

C 

1 .830676-01 

3.999446-01 

11 

G 

-1 .677486-01 

-3.770496-01 

12 

G 

1 .677486-01 

3 .770496-01 

13 

G 

-2.066426-01 

-4.734196-01 

14 

G 

2.066426-01 

4.734196-01 

15 

G 

-1 .841036-01 

-4  244466-01 

16 

G 

1 .841036-01 

4.244466-01 

17 

G 

-1 .645956-01 

-3.867446-01 

16 

G 

1 .64595E-01 

3.86744E-01 

19 

G 

-1 .491066-01 

-3.213946-01 

20 

G 

1 .491086-01 

3.213946-01 

21 

G 

-1 .66272E-01 

-3.787396-01 

22 

C 

1 .662726-01 

3.787396-01 

23 

G 

-1 .926596-01 

-4.650746-01 

24 

G 

1 .92659E-01 

4.660746-01 

25 

G 

-1 .67187E-C1 

-4.113766-01 

26 

0 

1 .67187E-01 

4.113766-01 

27 

G 

-1 ,335476-01 

-3. 858166-01 

28 

G 

1  33547E-01 

3.858166-01 

29 

G 

-1 .438916-01 

-2.96662E-01 

30 

C 

1.438916-01 

2.96662E-01 

31 

G 

-1 .614526-01 

-3.446206-01 

32 

G 

1 .514626-01 

3.446206-01 

33 

G 

•1  .676516-01 

-4.259786-01 

34 

G 

1  .675616-01 

4.259786-01 

36 

G 

-1.311166-01 

-3.761616-01 

36 

G 

1 .311166-01 

3.761616-01 

37 

G 

-1 .064616-01 

-3.666686-01 

38 

G 

1 .084616-01 

3.666686-01 

39 

G 

-1.274746-01 

-2.494206-01 

40 

G 

1.274746-01 

2.494206-01 

41 

G 

-1 .302S3E-01 

-2.918766-01 

42 

G 

1.302536-01 

2.918766-01 

43 

G 

-1.392376-01 

-3.653736-01 

44 

G 

1.392376-01 

3.653736-01 

4S 

G 

-9. 864646-02 

-3,171226-01 

46 

G 

9.864646-02 

3.171226-01 

47 

0 

-7.501166-02 

-3.044566-01 

48 

G 

7.501156-02 

3.044666-01 

49 

G 

-1.056266-01 

-1.801626-01 

50 

G 

1 .066266-01 

1.061526-01 

61 

G 

-1.037316-01 

-2.100966-01 

62 

G 

1 .087316-01 

2.190966-01 

83 

G 

•1. 080066^1 

-2 . 706976-01 

64 

Q 

1.010066-01 

t , 700076-01 

68 

G 

-0.230096-02 

-2.127006-01 

56 

G 

6.238086-02 

2.127006-01 

57 

G 

-4.138696-02 

-2.324926-01 

58 

G 

4.138696-02 

2.324926-01 

59 

G 

-7.970666-02 

-1.161016-01 

60 

G 

7.970666-02 

1.161016-01 

61 

G 

-7.472666-02 

-1.392416-01 

62 

G 

7.472666-02 

1.392416-01 

63 

G 

-7 .063926-02 

-1 .742686-01 

64 

G 

7.063926-02 

1.742686-01 

65 

G 

-3.280776-02 

-1.326686-01 

66 

G 

3.280776-02 

1.326686-01 

67 

G 

-1.421936-02 

-1.427216-01 

68 

G 

1.421936-02 

1 .427216-01 

69 

G 

-4.112666-02 

-4.112766-02 

70 

G 

4.112866-02 

4.112766-02 

71 

-4.489166-02 

-6.971186-02 

72 

G 

4.489166-02 

6.971186-02 

73 

G 

-4.464606-02 

-8.836226-02 

74 

G 

4.464606-02 

8.836226-02 

75 

G 

-1 .932936-02 

-6.866036-02 

76 

G 

1.932936-02 

6.856036-02 

77 

G 

-6.268656-03 

-8.077946-02 

78 

G 

6.266656-03 

8.077946-02 

79 

G 

0.000006*00 

0.000006*00 

80 

G 

0.000006*00 

0.000006*00 

81 

G 

0.000006*00 

0.000006*00 

82 

G 

0.000006*00 

0.000006*00 

83 

G 

0.000006*00 

0.000006*00 

84 

G 

0.000006*00 

0.000006*00 

bL. 

G 

0.000006*00 

0.000006*00 

86 

G 

0.000006*00 

0.000006*00 

87 

G 

0.000006*00 

0 .000006*00 

88 

0 

0.0OOO0E*00 

0.000006*00 

ASTROS  VERSION  2  6/29/86 

STATICS  ANALYSIS:  BOUNDARY  1,  SUBCASE  1 


ENT  VECTOR 

T3  Rl  R2  R3 

1.60727E«01  0.e0000E‘»00  e.e000eE«00  0.00000E«00 

l.S0727E«01  0.00000E«00  e.e0000£«00  0.00000E«00 

i.6050iE«0i  0.00e00E«e0  0.e00e0£*e0  0.ee00eE«00 

i.6060iE«0i  e.0ee00E«00  0.ee000E<*00  0.e0000E«00 

1.704B0E«01  0.0ee00E-»00  0.e000eE«00  0.00000E«00 

1.70466£«01  0.00000E*00  0.00000E«00  0.00000E*00 

l.ei036E«01  0.000e0E«00  0.e0000E«00  0.00000E«00 

l.ei036E«01  0.00000E«00  0.ee0e0E*00  0.00000E*00 

1.92361E«01  e.00000E*00  0.00000E*00  0.00000E«00 

1.923S1E«01  0.00000E«00  0.00000E«00  0.e0000E«00 

i.6n8e£«0i  e.00ee0E*00  0.e0000E*00  0.00000E«00 

1.61186E«01  0.00000£«00  0.00000E«00  0.00000£«00 

1.61485E*01  0.00000E«ee  0.00000E«00  0.00000E«00 

l.B14e6E«01  0.00e00E«00  0.00000E*00  0.00000E«00 

1.61876E«01  0.00000E«00  0.00000E«00  0.00000E«00 

].6ie7SE«01  0.00000E*0e  0.0e000E*00  0.00000E«00 

1.61960E«01  0.e0000E«00  0.0e000E«00  0.00000E*00 

1.61950E«01  0.00000E«00  0.e0000E*00  0.00000E«00 

1.08286E«01  0.00000E«00  0.00000E«00  0.00000E«00 

1  0e286E*01  0.00000E»00  0.00000E«00  0.00000E«00 

1.0892BE*01  0.00000E*00  0.e0000E<»00  0.00000E«00 

i.ee92BE«0i  0.0e000E«00  0.ee00eE<»00  0.0e000E-»00 

l.e9003E«01  0.00000E«e0  e.00000E«00  0.00000E*00 

1.09003E«01  0.00000e«00  0.e0000E*00  0.00000E*00 

1.09049E«01  0.00000E«00  0.00000E*00  0.00000E*00 

1.09049E«01  0.00000E«00  0.000e0£«00  0.00000E*00 

].0823)£«01  0.00000E«00  0.00000E«00  0.0000eE«00 

1.0e23lE40I  0.00000E«00  0.00000E«00  0.00000E«00 

7.282B7E«00  0.00000E-»00  0.00000E*00  0.00000E«00 

7.282B7E«00  0.00000E*00  0.00000E«00  0.00000e«00 

7.34516E«00  0.00000E*00  0.00000E«00  0.00000E»00 

7.34B16E-»00  0.O0000E«00  0.00000E«00  0.00000E4^00 

7.32498E«00  0.00000E«00  0.00000E«O0  0.00000E«00 

7.32498E«00  0.00000E«00  0.0e000E«00  0.O0000E«00 

7.27819E*00  0.0e000E«e0  0.000e0E«00  0.O0000E*00 

7.27ei9E«00  0.00000E<»00  O.00000E«00  0.00000E«00 

7.i2e04E400  0.00000E40e  0.ee000E«00  0.00e00E«00 

7.12a04E«00  0.0e000E«00  0.000e0E«00  0.00e00E*00 

4.429e9E400  O.00000E«e0  0.e0000E400  0.e0000E«00 

4.42989E400  0.e0000E-»00  0.00000E«00  0.00000E400 

4.49267E400  e.00000E*00  0.00000E400  0.00000E«00 

4.49267£«00  e.00000E«00  0.e0000E400  0.00000E400 

4.4620eE«00  0.00000E^00  e.00000E*00  0.00000E400 

4. 402006400  0.00000E400  0.00000E400  0.00000E400 

4.184B9E400  0.00000E400  0.000006400  0.00000E*00 

4.304B9E400  0.00000E400  0.00000E400  0.00000E400 

4.17676E400  0.00000E400  0.00000E400  0.000006400 

4.17676E400  0.00000E400  0.00000E400  0.00000E400 

2.32284E400  0.00000E400  0.00000E400  0.00000E400 

2.82284E400  0.00000E400  0.00000E400  0.00000E400 

2.87864E400  0.00000C400  0.00000E400  0.00000C400 

2.37304E400  0.00000C400  0.00000C400  0.00000E400 

2.88001E400  0.00000E400  0.B0000C  00  0. 000006400 

2.80001E400  0.00000C400  0.00000E400  0.00000C400 

2.240076*00  0.000006*00  0.000006*00  8.000006*00 

2.240976*00  0.000006*00  0.00000E400  0. 000006*00 

2.00472E«00  0.000006*00  0.000006*00  0.000006*00 

2.00472E400  0.000006*00  0.000006*00  0.000006*00 

9.393966>01  0.000006*00  0.000006*00  0.000006*00 

9.39396E>01  0.000006*00  0.000006*00  0.000006*00 

9.70B186-01  0.000006*00  0.000006*00  0.000006*00 

9.70B16E>01  0.000006*00  0.000006*00  0.000006*00 

9.263476-01  0.000006*00  0.000006*00  0.000006*00 

9.263476-01  0.000006*00  0.000006*00  0.000006*00 

8.333686-01  0.000006*00  0.000006*00  0.000006*00 

8.333686-01  0.000006*00  0.000006*00  0.000006*00 

6.980046-01  0.000006*00  0.000006*00  0.000006*00 

6.980046-01  0.000006*00  0.000006*00  0.000006*00 

1.609436-01  0.000006*00  0.000006*00  0.000006*00 

1.609436-01  0.000006*00  0.000006*00  0.000006*00 

2.180816-01  0.000006*00  0.000006*00  0.000006*00 

2.180816-01  0.000006*00  0.000006*00  0.000006*00 

2.683686-01  0.000006*00  0.000006*00  0.000006*00 

2.683686-01  0.000006*00  0.000006*00  0.000006*00 

2.816716-01  0.000006*00  0.000006*00  0.000006*00 

2.81671E-01  0.000006*00  0.000006*00  0.000006*00 

1.876896-01  0.000006*00  0.000006*00  0.000006*00 

1.876896-01  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 

0.000006*00  0.000006*00  0.000006*00  0.000006*00 
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STATICS  ANALYSIS:  BOUNDARY  I,  SUBCASE  1 

STRESSES  IN  QUADRILATERAL  MEMBRANES  (QDMEMl) 

ELEMENT  LA^ER  STRESSES  IN  ELEMENT  COORD  SYSTEM  PRINCIPAL  PRINCIPAL  STRESSES  MAX 

ID  NO  NORMAL-X  NORMAL-Y  SHEAR-XY  STRESS  ANGLE  MAJOR  MINOR  SHEAR 


3  0  -2.493564E-»04  -6 . 225000E^0J  2.612795E404 

4  0  2.493564E*04  S.22&000E*03  -2  61279BE«04 

?  0-1  391456E*04  -3.417625E*04  8.3S6375E*03 

6  0  1 . 39)456E«04  3.417526E*04  -8 . 356375E*03 

7  0  3.43840eE*04  -2.278972E*04  3.294600£*04 

8  0  3.438406E*04  2.278972E404  -3 . 294600E«04 

9  0  -7 .008719E*03  -4 . 979883E*04  -  1 . 288219E *04 

10  0  7.008719e*03  4  979883E*04  1.280219E*04 

11  0  -0  916969E*03  -6 , 322420E-*04  1.304433E*04 

12  0  0.915969E*03  S.322420E«e4  - 1 . 304433E*04 

13  0  -9 . 701406E«03  -6 . 370888E«04  1.399e76E«03 

14  0  9  701406E*03  6.370e8eE«04  - 1 . 399876E*03 

15  0  -6. 123603E*03  -6 . 9401 28E-*04  1  .630808E-*04 

16  0  6  123563E*03  5.948128E*04  >1 .630808E«04 

17  0  -2 .097320E*03  -5 . 881309e*04  - 7 . 810977£*03 

18  0  2.097320E*03  6.881309E*04  7  810977E*03 

19  0  -5.921461E*03  -  7 . 21423 1 E *04  8 , 163492E*03 

20  0  5.921461E*03  7.21423lE*04  -0 . 163492E*03 

21  0  -4 ,059875E*03  -  7 . 454 122E*04  -3 .056547E*03 

22  0  4 ,059075E-»03  7 . 454 122E-*04  3  066547E*03 

23  0  -I ,0336O9E-03  -6 . 862981E*04  1.379294E*04 

24  0  1,033609E*03  6.862981E404  -1 . 379294E404 

25  0  -1 .9356645*03  -6 .905132E404  -8 . 192309E403 

26  0  1.935664E*03  6.906132E*04  8.192309E403 

27  0  -2 .925473E*03  - 7 . 907642e*04  6.127383E403 

28  0  2.925473E*03  7.907642E*04  -6 . 127303E-03 

29  0  -2 . 729094E*03  -8  336064E*04  -6.8946396*03 

30  0  2.729094E*03  8.3360645*04  6  894S39E*03 

31  0  -1.0149066*03  -7.9694146*04  1.192S90E*04 

32  0  1,014906E*03  7.9694146*04  -1.1925986*04 

33  0  -3  6203916*03  -7.6665546*04  -6  7072506*03 

34  0  3.6283916*03  7.6566546*04  6.7072606*03 

35  0  -5  0542426*03  -9  0330526*04  3  6912736*03 

36  0  5.0542426*03  9.0330526*04  -3 .691273E*03 

37  0  -3.2616086*03  -9.6620566*04  -9.8323286*03 

38  0  3.2615086*03  9.6620566*04  9.8323286*03 

39  0  -6  8483206*02  -8.8447676*04  9.0742866*03 

40  0  6.8483206*02  8.8447676*04  -9.0742856*03 

41  0  -5.3239326*03  -7,6640776*04  -4.7303896*03 

42  0  5.3239326*03  7.6640776*04  4.730389E-03 

43  0  -8  316775E«03  -9.1363886*04  -7.9666026*01 

44  0  8.3167766*03  9  1363886*04  7.9666026*01 

45  0  -4  8685636*03  -9.6693296*04  -1.4761686*04 

46  0  4.8685636*03  9.6693296*04  1.4761886*04 

47  0  -47618556*02-9.0746136*04  4.8164226*03 

48  0  4  7618666*02  9.0745136*04  -4 . 8l6422E«e3 

49  0  -1.2161786*04  -7.107522e-*04  -1.9694  716*03 

50  0  1.2161786*04  7.1076226*04  1.9694716*03 

51  0  -1.4774456*04  -8.7031276*04  -2.8797696*03 

62  0  1.4774466*04  8.7031276*04  2. 8797696*03 

53  0  -1.1633066*04  -9.6429276*04  -1.9994006*04 

64  0  1.1633066*04  9.6429276*04  1.9994006*04 

56  0  -6.8223936*03  -1.0070226*06  -3.7607096*03 

63  0  6.8223936*03  1.0078226*06  3.7687096*03 

67  0  -1.1623016*04  -4.0218026*04  1.0201696*04 

60  0  1.1623016*04  4.8218026*04  -1.0201096*04 

69  0  -1.9217326*04  -8.1039736*04  1.0634696*04 

60  0  1.9217326*04  8.1039736*04  -1.0034696*04 

61  0  -1,6447446*04  -9.3681796*04  2.1690736*03 

62  0  1.6447446*04  9.3601796*04-2.1690736*03 

63  0  -2.2106696*04  -1.0826016*05  3.1968766*03 

64  0  2.2106696*04  1.0826816*05  -3.1968766*03 


66.3331  1.2644636*04  -4.3005176*04  2.7924866*04 

-34.6689  4.3006176*04  -1.2844636*04  27924866*04 

19.7594  -1.0912786*04  -3.7177046*04  1.313213E*04 

-70.2406  3.7177046*04  1.0912786*04  1.3132136*04 

49.9898  4.0652686*03  -6.2039046*04  3.3452156*04 

-40.0102  6.2039046*04  -4.8662686*03  3.3462166*04 

-16.6283  -3.4298066*03  -6.3377746*04  2.4973976*04 

74.4737  6.3377746*04  3.4298066*03  2.4973976*04 

15  2448  -6.3609566*03  -6.6779226*04  2  570913E*04 

-74.7662  6.6779226*04  6.3609566*03  2.6709136*04 

1.4838  -9.666145E*03  -6.3746146*04  2,7040006*04 

-88.6162  6.374614E*04  9.6661466*03  2-7040006*04 

15.7182  -1.6339886*03  -6.4070866*04  3.1268436*04 

-74.2818  6.4070866*04  1.5339886*03  3.1268436*04 

-7.8047  -1.8266896*03  -6.9883736*04  2. 9028526*04 

82.1963  6.9883736*04  1.8266896*03  2.902e5?E*04 

6.9170  -4.9323266*03  -7.3131456*04  3  4099666*04 

-83.0830  7.3131466*04  4.9323286*03  3.4099666*04 

-2.4777  -3.9276606*03  -7.4673446*04  3.6372896*04 

87.6223  7.4673446*04  3.9276606*03  3.6372896*04 

11.1001  1.6724926*03  -7.1336916*04  3.6604206*04 

-78.8999  7.1336916*04  -1.6724926*03  3.6504206*04 

-6.8596  -9.6016026*02  -7.0036836*04  3.4643346*04 

83.1406  7.0036836*04  9.6016026*02  3.4643346*04 

4.6711  -2.4366866*03  -7.9666316*04  3. 8664866*04 

-86.4289  7.9666316*04  2.4366866*03  3.8604866*04 

-4.1696  -2.3004026*03-8.3779346*04  4.0739476*04 

85.0404  8.3779346*04  2.3004026*03  4.0739476*04 

8.4325  7.5307036*02  -8.1462126*04  4.1107596*04 

-81.6676  8.1462126*04  -7.6307036*02  4.1107596*04 

-6.2113  -3.0166466*03  -7.7167296*04  3.7075326*04 

84.7887  7.7167296*04  3.0166466*03  3.7076326*04 

2.4739  -4.8947626*03  -9.0490016*04  4.2797636*04 

-87.6261  9.0490016*04  4.8947626*03  4.2797636*04 

-6.0162  -2.2262856*03  -9.6656776*04  4.7165756*04 

83.9838  9.6656776*04  2.2252866*03  4.7165756*04 

6.8418  2.4358696*02  -0.9376096*04  4.4809846*04 

-84.1682  6.9376096*04  -2.4368696*02  4.4809846*04 

-3.7783  -6.0116366*03  -7.6963166*04  3.6970826*04 

86.2217  7.6953166*04  6.0115366*03  3.6970826*04 

-0.0660  -8.3166966*03  -9.1363966*04  4.1523636*04 

89.9450  9.1363966*04  8.3166966*03  4.1523636*04 

-8.9062  -2.6668326*03  -9 . 9005026*04  4.8224106*04 

81.0938  9.9005026*04  2.5568326*03  4.8224106*04 

3.0449  -2.1903626*02  -9.1001286*04  4.5391136*04 

-86.9551  9.1001286*04  2.1903526*02  4.5391136*04 

-1.9122  -1.2086036*04  -7.1140986*04  2-9527486*04 

88  0878  7.1140986*04  1.2086036*04  2-9527486*04 

-2.2787  -1.4659866*04  -8.7146866*04  3.6243006*04 

87.7213  8.7146866*04  1.4659866*04  3.6243006*04 

-12.7421  -7.0117606*03  -9.9960686*04  4.6469416*04 

77.2679  9-9960686*04  7.0117606*03  4.6469416*04 

-2.2632  -6.0738406*03  -1.0093086*06  4.7828466*04 

§7.7368  1.0093086*05  6.8738406*03  4.7828466*04 

14.6706  -8.9712606*03  -6-0870386*04  2.0949666*04 

-76  4294  6.0870386*04  8-9712606*03  2. 0949566*04 

9.4928  -1.7439126*04  -8.2817936*84  3.2889416*04 

-90.5074  9,2817936*04  1.7439126*04  3.2689416*04 

1.6095  -1.6386496*04  -9.3642746*04  3.8628136*04 

-88.3905  9.3642746*04  1.6386496*04  3-8628136*04 

2.1222  -2.1988226*04  -1.0837666*06  4.3194186*04 

-87.8778  1.0837866*05  2.1988226*04  4.3194186*04 
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stresses  in  triangular  m 

ELEMENT  LAYER  STRESSES  IN  ELEMENT  COORD  SYSTEM  PRINCIPAL 

tD  NO.  NORMAL-X  NORMAL-Y  SHEAR-XY  STRESS  ANGLE 

1  0  1  323869E*04  -  1 . e20e7  6E<04  -8 . 4 leOSAE^aj  -17.8379 

2  0  -I . 323889E*04  1.020e7EE*04  8.4ie094E>03  72.1622 


STATICS  ANALYSIS;  BOIKJARY  i,  SUBCASE  1 
MBRANES  (TRmEU) 

PRINCIPAL  STRESSES  MAX 

MAJOR  MINOR  SHEAR 

1.694e93E»04  -1 .291B00E<04  1 . 443097E-04 

1.291B00E»04  -1 .B94693E«04  1.443097E»04 


INTERMEDIATE  COMPLEXITY  WING 


ASTROS  VERSION  2  6/29/88 


STATICS  ANALYSIS:  BOUWARY  I,  SUBCASE  1 


S 

TRESSES 

IN  SHEAR 

PANEL 

S  (SHEAR 

) 

ELEMENT 

MAX 

AVERAGE 

SAPETY 

ELEMENT 

MAX 

AVERAGE 

SAFETY 

10. 

SHEAR 

SHEAR 

MARGIN 

10. 

SHEAR 

SHEAR 

MARGIN 

6& 

9.210000E»03 

<7.9e6689E«03 

1 .0E430 

66 

6.624000E403 

-4.878268E403 

1 .0E430 

67 

I  .  710e92E«04 

-1  .5l0e96£«04 

1 .0E430 

68 

1  .e64156E*04 

'1.6076786404 

1 .06430 

69 

1  .  834200E«04 

1  .S64920£*04 

1 .0E«30 

70 

1 .963200E404 

1 .6917536404 

1  .06430 

71 

a  .  34e507E«03 

7  .553504E«03 

1 .0E430 

72 

2.237213E404 

1 .9893326404 

1  .06*30 

73 

4 .066000E«03 

3.449971£«03 

1 .0E430 

74 

1 .352300E404 

1.160491E«04 

1 .0E430 

7S 

1 . 2e73S0E«04 

>1  .  ie3700E«04 

1 .0E430 

76 

2.714266E403 

'2.409632E403 

1  .06430 

77 

8, 721S00E«03 

7.405316£«03 

1 .0E430 

78 

1.646800E404 

1 .4207B7E404 

1  .0E430 

79 

6.7e8029E«03 

-e.l31640E-8e3 

1 .0E430 

80 

4.7ie2nE403 

4.1838666403 

1.06430 

81 

1  .080700E-04 

9. 181279E-*03 

1 .0E430 

82 

1.749626E404 

1 .511384E404 

1 .0E430 

83 

2. 125704E*03 

-1  .gi9040E«03 

1 .0E430 

84 

1 .042047E404 

9.2387966403 

1  .0E430 

85 

2 . 330800E«04 

1 .981249E«04 

1 .0E430 

86 

2.612588E404 

2.2B8947E404 

1 .0E430 

87 

4 .084244E«03 

3.686403E«03 

1 .0E430 

88 

1  .851117E404 

1  .640405E404 

1 .0E430 

89 

3  848934E«04 

3.272966E«04 

1 .0E«30 

90 

3.676302E404 

3. 179983E404 

1 .06430 

91 

2. 189809E«04 

1 .974433E«04 

1 .0E«30 

92 

3.e57693E404 

3.4164346404 

1 .0E430 

93 

4 .024393E«04 

3.457490E«04 

1 .0E430 

94 

3.738680E404 

3.2687666404 

1 .06430 

95 

2.203178E-04 

1 . 9644536 «04 

1 .0E430 

96 

3.885516E404 

3.402794E404 

1 .06430 

97 

2.601600E-»04 

2.2ie626E«04 

1 .0E430 

96 

2.065000E404 

1 .655005E404 

1 .0E430 

99 

2.048700E«04 

1 .8S9643E«04 

1 .0E430 

100 

1  .e43281E404 

1  .6865666404 

1  .06430 

101 

1  .429676E*04 

1  .SieiTSE-eA 

1 .0E430 

102 

3.447203E403 

3.1654356403 

1 .0E430 

103 

1  .  436901E*04 

-1.32831l£*04 

1 .0E430 

104 

3.342700E404 

3.1663586404 

1 .0E430 

105 

3.299100E«04 

2.931614E«04 

1 .0E430 

106 

3.337d60E404 

3.0030136404 

1 .06430 

107 

3.385626E«04 

3.077031£«04 

1 .0E430 

108 

3.335009E404 

3.0563696404 

1 .06430 

109 

3.101364E404 

2.862389E*04 

1 .eE430 

110 

2.7269706404 

2.5624826404 

1.06430 

in 

2.621402E«04 

2.474693E«04 

1.0E430 

112 

2.0113008404 

1.8086226404 

1 .06430 

113 

3.751600E«04 

3.3420046404 

1 .0E430 

114 

3.674076E404 

3.3105166404 

1 .06*30 

115 

3.62e563E«04 

3.2999376404 

1 .0E430 

116 

3.633713E404 

3.241692E404 

1 .06*30 

117 

3.517e81E'*04 

3.2491366404 

1.0E430 

118 

3.728138E404 

3.6915126404 

1 .06*30 

119 

4 . 104376E«04 

3.9592786404 

1 .0E430 

INILKULDIATE  complexity  wing 


ASTROS  VERSION  2  6/29/88 


STATICS  ANALYSIS;  BOUtCARY  I,  SUBCASE  1 
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Appendix  D 

Random  Element  Selection 
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The  random  element  to  be  damaged  was  selected  thusly:  A 
uniform  random  number  generator  was  called  to  provide  a  number 
between  0  and  100.  The  seed  number  was  incremented  by  one 
every  time  the  batch  file  was  run.  If  the  number  was  between 
0  and  90,  a  quad  element  subroutine  was  called.  If  between  90 
and  99,  a  shear  panel  subroutine  was  called.  If  between  99  and 
100,  the  triangular  subroutine  was  called. 

The  quad  subroutine  called  another  random  number,  normalized 
to  0  to  31,  corresponding  to  the  number  of  even  quad  elements. 
This  number  was  multiplied  by  two,  to  give  an  even  number. 
Then,  that  number  had  2  added  to  it  to  compensate  for  the 
first  two  numbered  elements  not  being  quad  elements.  That 
number  and  that  number  minus  1,  were  considered  the  choice  of 
elements,  giving  an  even  numbered  quad  element  and  the  odd 
numbered  element  above  it.  These  numbers  were  returned  to  the 
section  of  the  program  which  determined  if  they  were  in  the 
section  of  the  wing  eligible  for  damage. 

The  shear  panel  subroutine  called  a  random  number,  normalized 
to  0  to  54,  corresponding  to  the  number  of  shear  panels.  This 
number  was  returned  to  the  section  of  the  program  which  deter¬ 
mined  if  it  was  in  a  section  eligible  for  damage. 

The  triangular  subroutine  merely  jumped  to  the  section  of  the 
program  which  checked  to  see  if  the  tip  was  the  eligible 
section . 
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Weight  Estimation  Calculations 


The  original  weight  was  estimated  by  assuming  structural 
material  possessed  a  unit  weight.  Thus,  the  cubic  volume  of 
the  structure  could  be  used  for  weight  calculations. 

Most  quad  elements  were  about  the  same  size,  and  approximately 
rectangular  in  shape.  Element  29,  from  the  middle  of  the 
wing,  was  selected  as  representative.  Using  its  nodal 
coordinates,  its  side  lengths  were  computed  to  be  about  8.4, 
13.8,  9.1,  and  13.8  units  long.  This  gave  the  representative 
element  an  area  of  about  120  square  units.  Then,  the  thick¬ 
nesses  of  all  quad  panels  were  summed  up  to  a  single  thickness 
and  multiplied  by  the  representative  area  to  give  about  480 
cubic  units. 

Next,  the  volume  of  the  shear  panels  was  estimated.  Since  the 
wing  had  a  regular  shape,  root,  mid  and  tip  shear  panels  could 
be  averaged  together  to  give  a  good  representative  panel 
length.  There  were  two  basic  panel  lengths.  Long  panel  ran 
from  the  root  to  the  tip;  short  panels  from  the  leading  to 
trailing  edge.  Average  long  panel  length  was  14  units;  aver¬ 
age  short  panels  were  9  units  long.  Average  panel  height  was 
4  units.  The  thicknesses  of  all  the  short  panels  wre  summed, 
then  multiplied  by  9  by  4  to  give  about  24  cubic  units. 
Similarly,  the  long  panels  summed  up  to  give  an  additional  70 
cubic  units.  Rod  elements  contributed,  due  to  their  slender¬ 
ness,  no  significant  volume  (0.05  c.u.)  In  total,  the  esti¬ 
mated  volume  of  the  wing's  structural  elements  v.'as  about  57  5 
cubic  units. 
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Then,  using  the  representative  quad  panel  for  area,  all  the 
strengthening  thicknesses  were  added  up,  minus  the  lightening 
thicknesses  (see  Figure  22),  to  give  an  additional  volume  of 
about  24  c.u.  Thus  the  increase  in  volume  (and  weight)  was 
about  4%. 
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Abstract 

Finite  element  analysis  has  been  used  as  a  design  tool  for  many 
years,  with  structural  reliability  being  ensured  through  use  of 
a  liberal  factor  of  safety.  Unfortunately,  the  safety  factor 
is  a  blanket  insurance  against  all  hazards,  and  a  designer  has 
no  way  to  optimize  a  structure  against  any  particular  hazard. 
This  is  particularly  troublesome  in  the  fields  of  aero/astro 
design,  where  every  bit  of  mass  must  serve  to  maximum  utility. 

The  method  of  Stochastic  Finite  Element  Analysis  allows  a 
designer  to  model  any  loading  or  hazard  condition  as  closely  to 
reality  as  desired  by  using  an  appropriate  probability  distri¬ 
bution  function.  Through  a  Monte  Carlo  simulation,  the  finite 
element  model  is  subjected  to  the  probability  functions.  The 
cumulative  output  is  analyzed  for  trends  in  failure  probability 
and  the  design  is  altered  to  enhance  its  reliability,  repeating 
the  process  until  the  desired  level  of  reliability  is  achieved. 
The  resulting  design  is  optimal  for  the  imposed  conditions,  and 
compared  to  a  structure  designed  with  a  traditional  factor  of 
safety  approach,  is  either  lighter  or  more  reliable.  This 
demonstration  revealed  that  for  similar  reliabilities,  a 
stochastically  designed  wing  was  20%  lighter  than  a  wing 
strengthened  by  the  factor  of  safety. 

The  major  drawback  in  applying  the  method  of  stochastic  finite 
element  analysis  is  that  very  large,  complex  models  can  require 
extraordinary  amounts  of  computer  resources. 
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